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^adulating  surface  and  a  surface  which  is  uniformly  rough,  i.e.  a 
pseudo  white  noise  surface,  are  considered  in  detail.)  An  approximate  ap¬ 
proach  is  presented  for  determining  the  mean  or  coherent  scattered  field 
produced  by  a  uniformly  rough  surface  and  the  resulting  solution  clearly 
shows  the  important  effects  of  multiple  scattering.  For  the  class  of 
surfaces  for  which  decorrelation  does  not  imply  statistical  independence, 
it  is  found  that  the  incoherent  power  scattered  by  such  surfaces  com¬ 
prises  two  parts.  The  first  part  is  the  conventional  diffuse  term  while 
the  second  part  is  specular  in  that  it  exists  only  at  the  specular 
scattering  angle.  For  a  gently  undulating  exponentially  distributed  sur¬ 
face,  the  incoherent  specular  power  is  found  to  exceed  the  coherent 
scattered  power  when  the  Rayleigh  parameter  is  larger  than  unity.  These 
results  have  significant  implications  with  regard  to  the  class  of  sur¬ 
face  statistics  to  be  considered  for  fitting  measured  surface  statistics. 
Finally,  the  far  field  approximation  as  it  applies  to  rough  surface 
scattering  is  examined  in  detail  for  arbitrarily  distributed  surfaces. 

The  failure  of  the  approximation  when  the  incident  illumination  is  un¬ 
bounded  is  explained.  For  a  bounded  incident  plane  wave,  the  far  field 
approximation  for  the  mean  scattered  field  is  shown  to  be  valid  provided 
that  the  illuminated  area  on  the  mean  surface  is  very  large  compared  to 
the  wavelength  squared.  The  validity  of  the  conventional  far  field  ap¬ 
proximation  for  the  variance  of  the  scattered  field  requires  a  large 
illuminated  area  which  also  encompasses  many  decorrelation  lengths  on  the 
surface.  The  results  obtained  relative  to  the  validity  of  the  far  field 
approximation  are  general  in  that  they  hold  for  any  degree  of  multiple 
scattering  on  the  surface. 
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1.0  INTRODUCTION 


The  overall  objective  of  this  research  is  to  provide  models  for  scatter¬ 
ing  from  randomly  rough  terrain  which  can  be  used  with  a  high  degree  of  con¬ 
fidence  in  systems  engineering  applications.  In  order  to  accomplish  this  goal, 
two  avenues  of  approach  are  actively  being  pursued.  The  first  approach  com¬ 
prises  the  refinement  of  existing  theories  while  the  second  approach  entails 
the  development  of  new  theories  which  go  beyond  existing  models. 

At  first  glance  it  might  appear  that  such  a  two-pronged  approach  does  not 
represent  an  efficient  use  of  resources;  however,  this  is  not  true.  For 
example,  one  may  argue  that  existing  theories  should  be  thoroughly  tested 
before  expending  efforts  on  advanced  techniques.  However,  the  very  nature  of 
the  problem  precludes  a  rigorous  and  complete  validation  of  existing  models 
because  these  models  are  based  upon  mathematical  simplifications  and  approxima¬ 
tion  which  are  not  fully  understood.  That  is,  a  full  understanding  of  the 
limitations  of  existing  models  will  only  come  about  when  these  models  have 
been  improved  upon  by  overcoming  some  of  their  fundamental  inadequacies.  Con¬ 
versely,  existing  models  should  be  developed  and  refined  to  the  fullest  extent 
possible  so  that  they  can  be  compared  to  measurements  for  the  purpose  of  moving 
toward  some  understanding  of  their  capabilities  and  limitations. 

In  research  of  this  nature  where  there  is  an  immediate  need  for  accurate 
models,  one  must  be  constantly  alert  to  the  pittfalls  of  using  inherently 
limited  models  simply  because  they  produce  numbers.  In  fact,  the  technical 
literature  on  rough  surface  scattering  is  permeated  with  the  application  of 
models  to  situations  which  are  totally  outside  the  range  of  validity  of  the 
models.  The  apparent  success  of  such  efforts  is  due  in  large  measure  to  the 
clever  selection  of  the  surface  statistics  rather  than  the  true  accuracy  of 
the  scattering  model.  This  has  led  to  a  great  deal  of  confusion  as  to  which 
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model  works  for  what  class  of  surfaces  and,  most  importantly,  why.  Since  it 
is  extremely  doubtful  that  the  kind  of  surface  measurements  necessary  as  inputs 
to  the  model  will  ever  be  available,  the  only  alternative  is  to  turn  to  theo¬ 
retical  approaches  for  improving  existing  models.  That  is,  for  scattering 
situations  where  there  is  significant  concern  that  the  existing  models  may  not 
be  valid,  it  would  appear  that  improvement  of  the  theory  is  the  only  viable 
approach.  For  this  reason,  the  majority  of  our  efforts  during  this  period  have 
been  devoted  toward  developing  and  improving  upon  our  fundamental  knowledge  of 
rough  surface  scattering. 

Very  often  when  one  conducts  research  directed  toward  fundamental  issues 
it  is  of  a  very  esoteric  nature  and  it  has  very  limited  practical  application. 
However,  such  is  not  the  case  with  this  research.  In  fact,  after  establishing 
some  very  basic  techniques,  it  is  shown  that  these  can  be  used  immediately  to 
derive  results  which  are  directly  applicable  to  real  problems.  For  example, 
one  often  sees  in  the  technical  literature  on  random  surface  scattering  use 
made  of  approximations  derived  from  deterministic  scattering  theory.  While 
such  approximations  must  hold  on  a  sample  by  sample  basis,  it  is  not  clear  that 
they  pass  unaltered  through  the  averaging  process  necessary  to  construct  the 
mean  scattered  field  and  power.  With  the  techniques  developed  herein,  it  is 
possible  to  follow  these  approximations  all  the  way  through  the  averaging  process 
and  therefore  rigorously  derive  their  impact  upon  the  statistical  moments  of 
the  scattered  field.  Thus,  while  some  of  the  topics  addressed  in  this  interim 
report  may  appear  to  be  far  removed  from  useful  terrain  scattering  models,  it 
must  be  remembered  that  they  are  the  beginnings  of  an  attempt  to  build  a  solid 
foundation  for  rough  surface  scattering  theory  which  has  heretofore  been  less 
than  satisfactory. 
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Summary  of  Results 

In  Section  2.0  an  integral  equation  for  the  stochastic  Fourier  transform 
of  the  current  induced  on  a  perfectly  conducting  rough  surface  by  an  incident 
plane  wave  is  developed.  The  stochastic  Fourier  transform  of  the  current  is 
defined  as  the  multidimensional  Fourier  transform  with  respect  to  all  stochastic 
surface  characteristics  upon  which  the  current  depends,  i.e.  it  is  a  transform 
from  Vn£,n-0,1,2,' • • ,  space  to  kn>n=0,l,2,’ • • ,  space  where  represents  the 
nch  order  directional  derivative  of  the  surface  height  £  .  The  integral  equa¬ 
tion  is  developed  by  multiplying  the  magnetic  field  integral  equation  for  the 
current  in  coordinate  space  by  a  Fourier  kernel  involving  all  orders  of  sur¬ 
face  height  derivatives  and  their  associated  transform  variables  and  then  aver¬ 
aging  this  equation.  By  converting  averages  over  V11^  to  convolutions  in  k^- 
space,  there  results  a  singular  integral  equation  of  the  first  kind  and  of 
infinite  dimensions.  The  merit  of  this  approach  centers  about  the  fact  that 
one  can  clearly  see,  in  the  integral  equation,  the  effect  of  all  the  higher 
order  surface  height  derivatives.  Furthermore,  it  is  possible  to  truncate 
the  infinite  dimensionality  of  the  integral  equation  by  retaining  only  those 
orders  of  surface  height  derivatives  which  have  a  significant  variance.  Rela¬ 
tionships  are  also  developed  which  show  how  the  stochastic  Fourier  transform 
of  the  current  may  be  employed  to  determine  the  statistical  moments  of  the 
scattered  field.  Limiting  cases  of  a  gently  undulating  surface  and  a  uniformly 
rough  or  pseudo  white  noise  surface  are  explored  relative  to  the  average  scattered 
field  generated  by  them.  For  the  uniformly  rough  surface,  an  exact  one  dimen¬ 
sional  singular  integral  equation  for  the  average  scattered  field  is  obtained 
and  it  is  found  to  be  very  similar  to  the  first  hierarchy  integral  equation 
resulting  from  the  diagrammatic  technique.  Contrary  to  perturbational  tech¬ 
niques,  this  result  does  not  predict  that  a  uniformly  rough  or  pseudo  white  noise 
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surface  will  act  as  a  perfect  reflector. 

In  Section  3.0  an  approximate  technique  is  developed  for  estimating  the 
average  field  scattered  by  a  uniformly  rough  surface.  The  analysis  is  based 
on  the  remarkable  similarity  between  the  dominant  terms  in  the  magnetic  field 
integral  equation  for  the  current  for  both  the  uniformly  rough  and  the  gently 
undulating  surfaces,  i.e.  the  only  essential  difference  relative  to  the  mean 
scattered  field  is  the  effect  of  the  grossly  different  surface  correlation 
lengths  for  the  two  cases.  This  similarity  suggests  that  if  the  small  corre¬ 
lation  length  associated  with  the  uniformly  rough  surface  could  be  mathematically 
introduced  into  the  description  of  the  scattering  associated  with  the  gently 
undulating  surface  and  if  the  latter  scattering  problem  could  be  solved  then 
it  should  be  a  good  approximation  to  the  uniformly  rough  surface  scattering 
process.  The  essence  of  the  approach  is  mathematical  similitude;  that  is,  the 
technique  in  itself  does  not  necessarily  have  any  physical  interpretation. 
Consequently,  it  is  absolutely  essential  that  the  results  are  capable  of  being 
put  into  one  to  one  correspondence  with  one's  physical  understanding  of  the 
scattering  process.  To  accomplish  the  introduction  of  the  artificial  corre¬ 
lation  length  into  the  scattering  description  for  the  gently  undulating  sur¬ 
face,  the  gently  undulating  surface  is  replaced  by  a  discrete  approximation 
comprising  large,  non-overlapping,  flat  areas  having  random  elevations  with 
respect  to  the  mean  planar  surface.  The  average  scattered  field  produced  by 
this  approximate  surface  is  found  to  depend  upon  the  degree  of  correlation 
between  adjacent  large  flat  areas.  It  is  at  this  point  in  the  development 
that  the  artificial  correlation  length  is  introduced,  i.e.  rather  than  taking 
the  large  areas  to  be  highly  correlated  corresponding  to  the  true  gently  undu¬ 
lating  surface  they  are,  instead,  taken  to  be  un correlated.  The  resulting 
average  scattered  field  is  subsequently  found  to  depend  upon  the  number  of 
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interacting  areas  and  their  height  variance.  Furthermore,  the  results  show  that 
the  number  of  interacting  areas  which  give  rise  to  a  maximum  scattered  field 
increases  with  surface  roughness.  Physically,  this  observation  corresponds  to 
the  fact  that  as  the  roughness  increases  so  does  t  e  degree  of  multiple  scat¬ 
tering  on  the  surface,  i.e.  an  increase  in  multiple  scattering  implies  an 
increase  in  the  number  of  effective  interacting  areas  on  the  discrete  artifi¬ 
cial  or  substitute  surface.  Thus,  the  results  of  the  analysis  do  indeed  appear 
to  have  a  very  good  physical  basis.  Comparison  of  numerical  results  with  ap¬ 
proximate  results  from  the  diagrammatic  approach  show  very  good  agreement  over 
the  range  of  Rayleigh  parameter  where  one  should  expect  agreement. 

In  Section  4.0  the  problem  of  scattering  from  surfaces  which  are  neither 
Gaussian  distributed  nor  do  they  become  statistically  independent  as  they 
become  decorrelated  is  considered.  This  is  a  very  important  practical  problem 
because  there  is  certainly  no  guarantee  that  real  terrain  is  adequately  repre¬ 
sented  by  a  surface  roughness  which  is  Gaussian  distributed.  The  stochastic 
Fourier  transform  of  the  current,  developed  in  Section  2.0,  is  used  to  derive 
exact  expressions  for  the  mean  and  variance  of  the  scattered  field.  For  sur¬ 
faces  in  which  decorrelation  does  not  imply  statistical  independence  it  is 
found  that  the  variance  of  the  scattered  field  or  the  incoherent  power  has  two 
distinct  parts.  The  first  part  is  the  conventional  so-called  diffuse  power 
which  is  related  to  the  spatial  Fourier  transform  of  the  two  point  joint  char¬ 
acteristic  function  for  the  random  surface  characteristics.  The  second  part 
of  the  incoherent  power  is  specular  in  its  angular  behavior  and  it  is  determined 
by  the  difference  between  the  two  point  joint  characteristic  functions  for  de- 
correlation  and  statistical  independence.  A  numerical  comparison  of  this 
specular  incoherent  power  term  with  the  mean  scattered  field  for  a  gently  undu¬ 
lating  exponentially  distributed  surface  is  also  presented.  When  the  Rayleigh 
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parameter  is  much  smaller  than  unity,  the  power  associated  with  the  mean  scat¬ 
tered  field  is  the  larger  quantity.  Near  a  Rayleigh  parameter  value  of  one, 
the  two  quantities  are  nearly  equal.  When  the  Rayleigh  parameter  greatly 
exceeds  one,  the  incoherent  specular  power  becomes  much  larger  than  the  power 
associated  with  the  mean  scattered  field.  The  importance  of  these  results  to 
the  conventional  interpretation  of  scattering  measurements  are  discussed  in 
depth.  These  results  also  indicate  that  extreme  caution  should  be  used  in 
fitting  measured  surface  statistics  to  functional  forms  for  the  probability 
density  function  in  which  decorrelation  does  not  imply  statistical  indepen¬ 
dence  . 

Section  5.0  considers  the  applicability  of  the  far  field  approximation  to 
scattering  from  infinite  rough  surfaces  when  the  incident  plane  wave  illumina¬ 
tion  is  either  bounded  or  unbounded.  This  also  is  a  very  important  theoreti¬ 
cal  and  practical  problem  because  rough  surface  scattering  measurements  are 
always  interpreted  using  this  approximation  and  yet  its  validity  has  only  been 
partially  demonstrated  for  very  special  surfaces.  The  analysis  once  again 
uses  the  stochastic  Fourier  transform  of  the  current  in  conjunction  with  the 
exact  expression  for  the  scattered  field.  For  the  case  of  an  incident  plane 
wave  which  is  unbounded,  it  Is  shown  that  the  mean  or  average  scattered  field 
is  a  redirected  plane  wave  propagating  in  the  specular  direction.  The  ampli¬ 
tude  of  this  plane  wave  is  attenuated  relative  to  the  incident  field  by  the 
effects  of  the  surface  roughness.  Furthermore,  its  polarization  is  dependent 
upon  the  multiple  scattering  processes  occurring  on  the  surface.  The  variance 
of  the  scattered  field  or  the  incoherent  power  in  any  one  direction  is  found 
to  depend  upon  a  weighted  average  of  the  angular  spectrum  over  all  angles  (both 
visible  and  invisible)  or  directions.  This  result  clearly  shows  that  the  con¬ 
cept  of  a  far  field  is  not  applicable  to  the  case  of  unbounded  illumination 
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because  in  the  far  field  approximation  the  scattered  power  should  be  d  ter- 
mined  by  one  and  only  one  Fourier  component  of  the  angular  spectrum.  When 
the  incident  illumination  is  bounded,  the  far  field  approximation  for  the  mean 
scattered  field  is  found  to  be  valid  provided  the  cross  sectional  area  of  the 
incident  field  is  very  large  compared  to  the  ^cuare  of  the  electromagnetic 
wavelength.  This  caveat  is  necessary  to  insure  that  the  difference  between 
the  support  of  the  incident  illumination  and  the  support  of  the  current  is  not 
significant.  These  two  supports  will  be  different  when  multiple  scattering  on 
the  surface  is  important.  For  the  variance  of  the  scattered  field,  it  is 
further  necessary  to  assume  that  the  support  of  the  incident  illuminaiton  encom¬ 
passes  many  surface  decorrelation  intervals  in  order  to  bring  the  exact  result 
into  agreement  with  the  far  field  approximation.  Finally,  the  exact  results 
are  used  to  show  why  the  far  field  approximation  breaks  down  when  the  incident 
illumination  is  unbounded. 

All  of  the  following  material  has  either  been  accepted  for  publication  in 
journals  or  it  is  presently  in  the  review  process  for  publication  in  a  journal. 
Consequently,  the  style  of  the  sections  is  tailored  to  the  demands  of  the  par¬ 
ticular  journal.  In  order  to  avoid  a  complete  retype  of  the  material,  it  was 
decided  to  include  it  in  this  report  as  it  was  sent  to  a  journal.  This  choice, 
hopefully,  will  not  be  too  confusing  to  the  reader. 
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Abstract 

An  exact  alternative  approach  to  the  diagrammatic  technique  for  treating 
scattering  from  rough  surfaces  is  developed.  The  magnetic  field  integral  equa¬ 
tion  for  the  current  induced  on  the  rough  perfectly  conducting  surface  is  multi¬ 
plied  by  a  Fourier  kernel  involving  all  orders  of  surface  height  derivatives  and 
their  associated  transform  variables.  Averages  of  this  weighted  equation  are 
converted  to  convolutions  in  the  transform  domain.  The  result  of  this  opera¬ 
tion  is  a  singular  integral  equation  of  the  first  kind  of  infinite  diirensions 
(because  of  the  infinite  number  of  height  derivatives)  for  the  stochastic  Fourier 
transform  of  the  current.  A  procedure  is  developed  for  estimating  the  effects 

of  ignoring  one  or  more  surface  height  derivatives  in  terms  of  the  range  of 
validity  of  the  resulting  approximate  solution.  Special  limiting  cases  of 
very  gently  undulating  surfaces  and  uniformly  rough  surfaces  are  examined. 

New  and  illuminating  results  are  obtained  for  the  latter  case. 

Introduction  and  Summary 

There  is  a  very  definite  need  to  move  beyond  the  approximate  rough  surface 
scattering  solutions  provided  by  physical  optics  and  boundary  perturbation  the¬ 
ory.  The  first  clear  indication  of  this  need  was  the  failure  of  physical  optics 
and  perturbation  theory  to  accurately  predict  the  magnitude  of  the  average  field 
[1]  and  acoustic  pressure  [2]  scattered  by  a  roughened  water  surface  for  large 
surface  roughness.  Measurements  of  the  incoherent  average  scattered  power  [3] 
have  also  provided  results  which  are  not  in  agreement  with  an  appropriate 
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combination  of  the  above  noted  single  scattering  approximations.  Of  course  it 
is  one  thing  to  recognize  the  need  for  a  less  approximate  solution  but  quite 
another  to  obtain  such  a  solution! 

Recently,  a  significant  advance  toward  this  goal  has  been  provided  by  the 
work  of  Zipfel  and  DeSanto  [4]  and  subsequent  work  by  DeSanto  [5,6].  The  re¬ 
sults  obtained  in  these  analyses  were  based  upon  two  very  important  points 
which  should  be  emphasized.  First,  by  recasting  the  coordinate  space  integral 
equation  in  the  transform  domain,  the  authors  were  able  to  circumvent  problems 
associated  with  the  stochastic  nature  of  the  coordinate  space  integral  equation 
for  the  "surface"  Green's  function.  Second,  by  using  cluster  decomposition 
and  partial  summation  techniques  it  was  possible  to  express  the  solution  of 
the  transform  space  integral  equation  in  terms  of  integral  equations  involving 
successively  higher  orders  of  surface  and  field  interaction.  Numerical  solu¬ 
tion  of  the  integral  equation  involving  the  first  order  interaction  for  the 

average  scattered  field  (the  zeroth  order  interaction  being  the  physical  optics 
approximation)  has  clearly  shown  the  physical  optics  approximation  to  be  inade¬ 
quate  for  large  surface  roughness  and  small  correlation  length  [7].  The  strength 
of  this  result  lay  primarily  in  the  rigorous  nature  of  the  approach  and,  con¬ 
sequently,  the  generality  of  the  solution.  Unfortunately,  the  solution  does 
not  lend  itself  to  an  obvious  physical  interpretation.  Furthermore,  it  is  also 
not  always  easy  co  see  the  direct  effect  of  the  surface's  statistical  charac¬ 
teristics  in  the  higher  order  integral  equations. 

The  purpose  of  this  paper  is  to  present  an  approach  and  a  solution  which 
provides  a  bit  more  physical  insight  into  certain  aspects  of  the  problem.  In 
particular,  one  of  the  strong  points  of  this  approach  is  that  it  allows  a  direct 
correlation  between  surface  statistical  approximations  and  the  range  of  validity 
of  the  resulting  solution.  This  is  an  important  result  because  it  allows  one 


to  estimate  when  an  approximate  solution  is  valid  and  when  it  breaks  down  based 


upon  the  statistical  characteristics  of  the  rough  surface. 

The  analysis  is  based  upon  the  magnetic  field  integral  equation  (MFIE) 

for  the  current  induced  on  the  surface  of  a  perfectly  conducting  body.  Since 

any  statistical  moment  of  the  scattered  field  can  be  obtained  from  knowledge 

of  the  current,  this  is  a  logical  starting  point.  A  key  element  in  this  or  any 

analysis  Involving  stochastic  integral  equations  is  the  technique  employed  to 

overcome  the  stochastic  nature  of  the  integral  equation.  The  approach  used 

here  is  to  multiply  the  MFIE  by  the  Fourier  kernel  expfjk^C  +  j  I  k^+^V11?;} 

n=l 

(where  C(x,y)  is  the  stochastic  surface  height,  V£(x,y)  are  the  stochastic 
surface  slopes,  etc.)  and  then  average  the  resulting  equation  over  £,  V£,  etc. 
The  averages,  however,  are  converted  to  convolutions  of  the  terms  in  the  inte¬ 
gral  equation  and  the  Fourier  transform  of  the  appropriate  joint  probability 
density  function  for  the  stochastic  variables  £,  V£,  etc.  Thus,  whereas  Zipfel 
and  DeSanto  [4]  used  a  coordinate  space  transform  to  overcome  the  difficulties 
associated  with  the  stochastic  nature  of  the  problem,  this  approach  transforms 
all  of  the  stochastic  variables  (5,  V£,  etc.)  into  a  non-stochastic  space 
(k^,  k^,  etc.).  The  formal  result  of  these  and  other  manipulations  are  coupled 
integral  equations  for  the  x  and  y  components  of  the  current  which  are  Fredholm 
equations  of  the  first  kind.  The  dimensionality  of  these  two  integral  equations 
is  formally  infinite  because  the  current  depends,  in  general,  upon  C,  V£,  etc. 
The  explicit  dependence  of  the  integral  equations  upon  the  surface  statistics 
is  contained  in  the  single  and  two-point  joint  characteristic  functions  for 

the  variables  C,  V£,  etc.  By  examining  the  way  in  which  the  surface  statistics 
2  2 

(<C  >»  etc.)  appear  in  the  characteristic  functions,  it  is  possible  to 

translate  the  neglect  of  higher  order  surface  height  derivatives  directly  into 
criteria  for  the  range  of  validity  of  the  resulting  approximate  solution.  An 
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orderly  procedure  for  relating  mathematical  simplifications  to  the  surface 
statistics  is  presented. 

As  an  example  of  the  power  of  this  technique,  it  is  applied  to  surfaces 
having  a  correlation  length  much  larger  than  the  electromagnetic  wavelength 
(very  gently  undulating  surfaces)  and  to  surfaces  having  a  correlation  length 
much  smaller  than  a  wavelength  (uniformly  rough  surfaces) .  For  very  gently 
undulating  surfaces,  the  technique  provides  results  which  are  in  agreement  with 
the  scalar  single  scattering  or  physical  optics  approximation.  For  uniformly 
rough  surfaces,  the  stochastic  Fourier  transform  approach  is  found  to  yield 
two  parts  to  the  slope  normalized  surface  current,  i.e.  Jg(r)/l  +  | VJ; [  where 
Jg  is  the  current.  One  part  is  found  to  be  a  function  of  only  the  stochastic 
surface  height  while  the  other  part  depends  on  the  surface  height  and  slopes 
but  in  such  a  manner  as  to  contribute  nothing  to  the  average  field  scattered 
by  the  surface.  An  integral  equation  for  the  average  scattered  field  is  derived 
for  the  uniformly  rough  surface  and  it  is  found  to  be  similar  to  the  approximate 
equation  obtained  previously  by  Zipfel  and  DeSanto  [4].  Finally,  it  is  argued 
that  the  gently  undulating  and  uniformly  rough  cases  provide  lower  and  upper 
bounds,  respectively,  on  the  average  scattered  field  for  a  fixed  height  variance. 
Approach  Rationale 

The  stochastic  surface  height  C(x,y)  is  defined  with  respect  to  the  z  =  0 
plane  and  is  taken  to  be  a  zero  mean  statistically  homogeneous 

process.  The  region  above  the  surface  is  free-space  while  the  medium  be¬ 
neath  the  interface  is  perfectly  conducting.  For  a  "sufficiently  well  behaved" 
surface,  the  current  3  induced  on  the  surface  by  an  incident  magnetic  field  H* 
is  given  by  the  so-called  magnetic  field  integral  equation  [8,  pg.  354]; 


((r)  -  2n(r)  x  H^r)  +  ~  J n(r)  x  [jg(ro)  x  VQg(|r  -  rj  )]  dSo  (1) 


where  n(r)  is  the  upward  directed  unit  normal  to  the  surface  and  g(|r-rQj)  is 


proportional  to  the  free  space  scalar  Green's  function,  i.e. 


‘  [-  V  - +  S]/^1  +  *  s2' 

8<|r-rQ|)  -  exp (-jko|r-ro|)  /  |r-rQ| 


and  k  -  2ir/X  is  the  free-space  wavenumber.  The  quantities  £  *  3£/3x  and 

O  O  X 

C  •  3£/3y  are  the  x  and  y  components  of  the  surface  slope.  Expanding  the 
double  cross  product  in  (1),  converting  the  surface  integration  to  an  integra- 

+  "  2 


1 


tion  over  the  z  *  0  plane  though  dS 
o  o 


2  2 

5  +  C  yields  the  following; 

x  y 


C  "  +  £  drt  ,  and  multiplying 
xo  yo  ° 


both  sides  of  (1)  by  vf; 

J(r)  -  2N(r)  *  H^(r)  +  i  f  |[5(r) •7o8)J(r<))  -  [!(?) •  J (?o> JV^g (2) 


where 


5<r)  -  v4  +  ?x2  +  cy2  k  \(r) 


N(r)  -  /  +  C  2"+  C2  ’  n(r) 
X  y 

dr^  -  dx  dy 
t  o  o 

o 


and  the  integration  is  over  the  entire  zq  m  0  plane.  Since  the  current  must 

/  "  V  “  of  *>■  -► 

be  tangential  to  the  surface  and  +  5^  >  0  ,  N(ro)*J(ro)  »  0  and  the 

z-component  of  the  current  may  be  expressed  in  terms  of  the  x  and  y  components 
as  follows; 


J  (r  )  ■  £  J(r)  +  £  J(r) 
z  o  x  x  o  y  y  o 


(3) 


Substituting  (3)  in  the  right  side  of  (2)  yields  two  coupled  integral  equations 
for  J^fr)  and  Jy(r).  The  coupling  is  a  consequence  of  the  term  (N(r) • J(rQ) ]VQg 
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which,  using  (3)  to  eliminate  J2^rQ)»  has  the  following  x  and  y-components; 


[N(r)-J(ro)] 


■=  - 


K'C, 

'  r 


)J  (r  ) 

x  o 


VS 


)J  (r  )\ 

,  y  °  f 


ls_ 


(4) 


where  qQ  is  xq  or  yQ  .  Thus,  it  is  the  x  and  y-components  of  (2)  along  with 
(3)  which  provide  the  essential  relationships  for  the  vector  components  of  the 
current . 

Having  found  the  current,  it  is  then  possible  to  determine  the  far-zone 

A 

scattered  electric  field  in  the  direction  k  and  at  a  distance  R_  from  the  z  =  0 

8  ° 

plane  via  the  following  relationship  [8,pg.  361] 


1 


k  n 
o  o 

4tt 


8<Vks  *  k. 


/ 


J (r)exp(j  k  *r)  dr 


(5) 


where  r  ”  xx +  yy  + £(x,y)z,  k  *k  x+k  y+k  z,  dr  =  dxdy  ,  and  the  inte- 

s  sx  sy  sz  t 

gration  is  over  the  z  =  0  plane.  Of  particular  importance  however  are  the 

■f  .  i  2 

statistical  moments  of  the  scattered  field  such  as  <E  >  and  <  E  >  where  <•> 

s  1  s 

denotes  an  average  over  all  random  quantities.  In  view  of  the  factor  exp(jkgzO 
in  (5),  it  would  do  no  good  to  compute  <J(r)>  from  (2)  and  (3)  because  the  aver¬ 
age  of  the  product  is  not,  in  general,  equal  to  the  product  of  the  averages. 

One  could  multiply  (2)  and  (3)  by  exp(j  kgzC)  and  average  the  resulting  equa- 

~y 

tions  but  the  result  would  only  be  useful  in  computing  <ES>  •  That  is,  higher 
order  field  moments  would  require  a  new  solution  of  some  appropriately  tailored 
form  of  (2)  and  (3).  This  type  of  reasoning  leads  to  the  conclusion  that  what 
is  needed  is  a  one-time  solution  of  (2)  and  (3)  which  could  be  used  in  (3)  or 
any  order  self-product  of  (S). 

The  average  of  (5)  entails  multiplying  by  the  probability  density  function 
(pdf)  for  all  stochastic  variables  upon  which  J  depends  and  then  integrating 
the  product  over  all  stochastic  variables,  e.g. 
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[(r)exp(j  kaz  0>  •  /'•/  •  * )exp (j  k^z  C)p^ 


•  dC  dVi;  •  •  •  (6) 

where  It  has  been  assumed  that  the  order  of  integration  [over  r  in  (5)]  and 
averaging  can  be  arbitrarily  interchanged.  The  surface  slopes  (V£) ,  curva¬ 
tures  (A),  rate  of  change  of  curvatures  (V^O ,  etc.  are  included  because,  in 
general,  the  current  is  a  function  of  these  stochastic  variables.  Starting 
first  with  the  integral  over  £  in  (6),  it  is  obvious  that  this  is  identical  to 
the  Fourier  transform  (with  respect  to  £)  of  the  product  of  the  current  and  the 
pdf.  Since  the  Fourier  transform  of  the  product  of  two  functions  is  the  con¬ 
volution  of  the  transforms,  the  C- integration  can  be  expressed  as  follows; 

CO 

c,...)p1(C,V?,...)exp(j  kgz  O  dC 


CO 

J  1(^t,k1,VC,V2c,...)P1(ksz-k1,V;,V2C,...)dk1  (7) 


where  the  tilde  denotes  the  Fourier  transform  with  respect  to  (  .  In  a  simi- 

2 

lar  manner,  the  integrations  over  the  slopes (Vc),  curvatures (V  0 ,  etc.  can 

n 

also  be  written  as  convolutions  so  that  (6)  becomes  (with  s  *  E  i  ) 

n  1=1 


<J(r)exp(j  k  C)>  *  lim  - 

sz  n-*»  (2tt  ) 


08”  /-/3<Vkl>£ 2 . 


•  P1(ksz-k1,-k2,...,-kn)dk1dk2...dkn  (8) 

where  k  is  an  n-dimensional  transform  vector  associated  with  Vn£  and  the  tilde 
n 

2 

implies  Fourier  transforms  with  respect  to  V  £,  etc.  If  the  second 


2 

moment  of  the  scattered  field,  <|Eg|  >  ,  is  desired,  it  may  be  obtained  from 

l 

knowledge  of  the  following  products ; 


<J  (?)J*(r')exp[Jk  -  lim 

9  9  82  n-H» 


J 

q 


(rfki 


(r'  ,-kJ  , •  •  •  ,-k'  p.(k  -k  ,-k  -k!,-k  ,-k' ,•  •*,-k  ,-k' ;r  -r' )dk.dk!  ••  *dk  dk' 
(  q  t  1  2  n  |  r2  sz  1  sz  1  2  2  n  n  t  t  11  n  n 


(9) 

where  *  denotes  the  complex  conjugate  and  P2  ( • )  is  the  two-point  joint  prob¬ 
ability  density  function.  In  obtaining  the  above  result,  use  has  been  made  of 
the  fact  that  F[  J  *(£,VC,  •  •  •)  ]  ■  ^^(-k^.-k^,- •  *)|*  where  F  denotes  the  Fourier 
transform  operation.  The  Fourier  transforms  of  the  single  and  two  point  joint 
probability  density  functions  are  the  single  and  two  point  joint  characteristic 
functions,  resnectively  [9,  pg.  254]. 

It  is  obvious  from  (8)  and  (9)  that  if  the  Fourier  transform  of  the  cur- 


*► 

rent  with  respect  to  the  stochastic  variables  can  be  found  then  <Eg>  and 

< } Eg |  >  may  be  determined  by  a  straightforward  integration  of  known  functions. 

Thus,  the  problem  reduces  to  solving  for  the  stochastic  Fourier  transform  of 

the  current  using  (2)  and  (3).  Unfortunately,  one  cannot  simply  take  the 

stochastic  Fourier  transform  of  both  sides  of  (2)  in  a  direct  manner;  such  an 

operation  must  be  accomplished  within  the  framework  of  probability  theory. 

That  is,  both  sides  of  (2)  are  first  multiplied  by  the  stochastic  Fourier  ker- 
00  _ 

nel  exp{jkxC  +  j  E  kn+l*'7^  and  then  averaged  over  all  stochastic  variables. 
n-1 

The  averaging  operation  weights  each  term  in  (2)  by  an  appropriate  probability 

density  function.  For  example,  the  left  side  of  (2)  and  the  source  term  on 

the  right  side  are  weighted  by  the  single  point  joint  pdf  while  the  integral 

term  is  multiplied  by  the  appropriate  two-point  joint  pdf.  The  averages  over 

C»VC,***  and  E  ,V  £  are  then  converted  to  integrations  with  respect  to 
0  0  0 
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the  transform  variables  by  means  of  convolution  identities.  The  net  result  is 
an  infinite  dimensional  Fredholm  integral  equation  of  the  first  kind  for  the 
stochastic  Fourier  transform  of  the  current.  This  is  the  essence  of  the  sto¬ 
chastic  Fourier  transform  technique. 

Detailed  Development 

While  the  basic  concept  of  the  stochastic  Fourier  transform  approach  is 
straightforward,  there  are  a  number  of  details  associated  with  converting  (2) 
into  the  appropriate  integral  equation  which  require  amplification.  Further¬ 
more,  there  is  at  least  one  intermediate  result  which  has  significant  implica¬ 
tions  and  should,  therefore,  be  derived  in  an  orderly  manner.  Since  the  con¬ 
version  of  (2)  requires  dealing  with  an  infinite  number  of  stochastic  variables 
(£,V£»V  £,*••)  and  their  associated  transform  coordinates  (k^,k2,k^,* * •) ,  the 
algebra  and  symbolish  is  very  tedious  especially  because  of  the  term 
[N(r)* J(ro)]VQg  .  Ignoring  this  term  in  (2)  does  not  impact  the  generality  of 

-►  -f  -V 

the  approach  because  if  one  understands  how  the  term  [N(r) J(r)  is  treated 
then  inclusion  of  the  neglected  term  can  be  accomplished  almost  by  inspection. 
Thus,  ignoring  the  coupling  term  in  (2)  yields  the  following  two  decoupled 
scalar  equations; 


Vr)  ■ 


2q*[N(r)  x 


n  (r)  ] 


*/{- 


I&,  _ 

3x 

o 


JL&.  +  2&A 

y  3y0  HJ  q 


J  (r  )dr 


(10) 


where  q  ■  x  or  y.  From  a  practical  point  of  view,  ignoring  the  coupling  term 

o 

in  (2)  implies  that  the  curvatures  (V  O  and  higher  order  derivatives  of  the 
surface  height,  according  to  (4),  are  negligibly  small;  however,  this  simplifi¬ 
cation  will  not  be  used  in  this  section. 

Multiplying  the  left  side  of  (10)  by  the  product  of  the  Fourier  kernel 
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} 


CO  ^ 

exp{jkxC  +  j  Z^  k^-V11  <;}  and  the  single  point  probability  density 
2 

Pi(S»^C»V  £,•••)  ,  integrating  over  all  the  stochastic  variables  to  form  the 
average,  and  converting  these  integrations  into  convolutions  in  the  Fourier 
transform  domain  yields 

<J  (?)  expire  +  J  r  [■■■[  WlV-'V 

n  (2tt)  n  J  J 


h  (kl  -ei*k2  *kn  '  *n}  d3ld^2  “  *  d^n 


(11) 


where  r  «  xx  +  yy  ,  J  (•)  is  the  stochastic  Fourier  transform  of  J  (r), 
c  q  q 

P^(’)  is  the  characteristic  function  for  the  single  point  joint  pdf,  i.e. 
the  stochastic  Fourier  transform  of  the  joint  pdf.  Taking  the  incident  mag¬ 
netic  field  on  the  surface  to  have  the  form 


H1^)  -  HQhexp(-jk1  *^t-jklzC) 


(12) 


the  term  2q*[N(r)  *H  (r)]  may  be  expressed  as  follows; 


2q‘[N(?)  xhV)]  -  2HQ[d+  Cjcx  +  cVlexp  (- jkt  -  j  ;) 


y  y 


(13) 


where  k^  ■  +  k^z  is  the  incident  wave  vector  of  magnitude  kQ  and  direc¬ 
tion  k^  .  The  factors  ,  and  are  determined  by  the  polarization  of 

the  Incident  magnetic  field,  i.e.  -  q*(z*h)  ,  *  -q*(xxfi)  , 

Z  X 

Cy  *  "q'(yx  h)  ,  and  are  independent  of  the  spatial  (rt)  and  stochastic 
variables.  The  average  of  the  product  of  (13)  and  the  Fourier 
kernel  is  a  straightforward  multidimensional  Fourier  transform,  i.e. 
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<2q*  [N(r)  x  H1^)}  exp  [jk  C  +  j  E  k  •7nC]> 

1  n«l  n+x 

-  2Hoe«p  (-jVvfc’  ♦  jc’3k^  + 


pl(kl  _kl*,k2»k3''  ”) 


where  k0  ■  k-  k„  +  k-  k_  .  The  derivatives  with  respect  to  k„  and  k_ 

2  2x  2x  2y  2y  2x  2y 

are  a  consequence  of  the  slope  terms  £  and  £  ,  respectively,  in  (13). 

x  y 

Averaging  the  integral  term  in  (10)  requires  the  general  two-point  prob¬ 
ability  density  function  P2(C»£0,7?,7  C0,7^£,V^C  ,...;r  -r  )  because  multi- 

o 

plication  by  the  Fourier  kernel  introduces  the  additional  stochastic  variables 
2  3 

7  £, 7  £,*••  which  are  not  contained  in  the  kernel  of  the  integral  equation 
(10).  Thus,  the  average  of  the  product  of  the  Fourier  kernel  and  the  integral 
term  in  (10)  may  be  written  as  follows; 


/"/(•  & '  Sr  r£  +  vVvv°''": 


‘*e  [iVi  +  i  JL  P2<c,c0,'^Vv"->d't  dt  dh m  dVoV 


Assuming  that  the  orders  of  integration  may  be  arbitrarily  interchanged,  the 
^-integration  can  be  written  as  a  convolution  of  the  ^-Fourier  transforms  of 
the  Green's  function  derivatives  and  P2(*)  •  Noting  that 


9g(rt-rt  ,C-C0> 


“  -  exPUW  FC 


3g(rt  “  rt  .O 

o 
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where  F  denotes  the  Fourier  transform  with  respect  to  ?  and  3q  is  the  trans¬ 
form  variable,  and  substituting 


g(Art,Bo)  *  (g(Art.O}  -  J  g(Art»C)  exp  (j3QC)  d? 


f  .O  \  f  3g(Ar  , 

gc(Art,&0)  -  F?  [  - ^ - J  =  J  - ^ 


C) 


exp(jB0C)  dC 


in  the  convolution  integrations  with  Ar  “  r  -  r  reduces  (15)  to  the  fol- 

o 

lowing  form; 


<•> 


"  /  '/  {' 


3g(Art,$o>  3g(Art,30) 


3x 


-  S  --if—  - 


f8o>} 


J  (r  ♦£  ,V  £  » *  •  •  )  exp 
q  t  o  o  o 


[jB 


0^0 


+  j 


CO 

E 

n=l 


n+1 


?nt] 


*  P2(kl‘Bo^o’V^*Vo^o'’";^t)  d6od^t  d?od7CdVoCo“'  (16) 

o 

where  the  tilde  symbol  denotes  the  Fourier  transform  of  P2(*)  with  respect 
to  £  . 

The  ^-integrations  in  (16)  may  be  expressed  as  convolutions  with  a  shift- 

2 

ed  argument  due  to  the  exponential  factor  in  (16).  The  ^QC0>  £0>-**  inte_ 

grations  may  be  represented  by  convolutions  with  no  shift  in  the  argument. 

2 

Finally,  the  integrations  over  V£,V  ?,•••  may  be  written  as  straightforward 
Fourier  transforms.  Accomplishing  all  of  these  operations  yields  the  follow¬ 
ing  result; 
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<•>  ■  lim 
n-*» 


(2tt) 


kf/( 


3g(Ar  ,B  )  3g(Ar  ,B  ) 

J— +  J  -ir~\  -s 

o  2x  'o  2y 


,<4VV} 


J  (rt  , ^ ,12 ,  •  •  •  »\)P2 (krS0 » Po~Bl,k2 *"^2 *  “  ’  ,kn’Ai ’ Art)drt  d3od^l '  ’ '  d^n 
o  o 

(17) 


where 


“  /'"/Jq(Vto,VoC,,'"’^o)eXP{j^0+j  j/n+l*VoCo} 


d £  dV  £  •  • .dV"^ 

O  O  O  0  0 


(18) 


and 


P2(kreo,eo“6l,k2,"^2’*”,kn’"^n;Art)  ~j  jP 2(5*VV^*Vo’‘ 


exp 


{ 


j(kx-Bo);  + 


J«W‘o 


+  j 


n 

£ 

1-1 


n+1 


■Vn£  - 


n 

£ 

1=1 


n+1 


■vnc 

0  O 


} 


•  d£  d£  dVC  dV  t,  •••dVn^dVnC  (19) 

0^00  o  o 

The  average  of  (10)  weighted  by  the  Fourier  kernel  is  thus  found  by  equating 

(11)  to  the  sum  of  (14)  and  (17).  However,  before  doing  this  the  substitu- 

**►  *►  “► 

tion  Ax  *=  x-x  and  Ay  *  y-y  or  Ar  =  r  -r  is  first  made  in  (17);  then  the 
o  o  ttt 

o 

equating  of  terms  yields 
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lim 
n-*»  (2tt) 


~%r  /■"/  3,<?fer52--"'6n)fi<kr6i-V52----V*„>d6id52---<,6„ 


2H 


“P  <-  J  *le-V[4  ♦  i  \2x  *  1  <$  \  J  ?!  <krtA  ’  V  ’ 


+  lim 


n-«°  (2tt) 


_  J  'f  f.  _«»vBo)  .  .  ,  agca?t,eo) 


rn  3k2x  *  J  sif"^  3k2y  + 


Jq(rt~Art,0l,02*‘'‘,en)  ei,k2,_B2’“  '  ,kn’“en:Art) 


•  dAr  dt$  d0,  •  •  -  dB 
t  o  1  n 


(20) 


Since  all  the  three  terms  in  (20)  must  exhibit  the  same  dependence  upon  r^  and 


only  J  and  exp(-jk  *r  )  are  functions  of  r  ,  (20)  implies  that  J  may 
written  in  the  following  factored  form; 


be 


‘Vrt,Sl’e2’"')  “  6Xp  ^  ki  *rt* 


(21) 


Substituting  this  result  in  (20)  and  rearranging  terms  produces  the  follow!.. £ 
integral  equation  for  (8^ ,62* *  *  * )  J 


lim 
n-H»  (2tt) 


^•n  J'"J  "  ,dn)  [ kl <kl'8l ,k2_B2 • "  ' 'k.rBn) 


-  r2(B1,i.l.62.k2>--.in.kI1)]  d6l  d|2---dln 


•) 


(22) 


where 
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*  -  i  f  ff  3l(Art,Bo)  3g(Art,6o) 

r2(61,kl’82,k2,'“  *\,kn)  “  (2it)2  J  J  3Ax  3k2x+J  3Ay  *k2y 


+  P2(kl'6o,So*Sl,k2,"^2’”  ’,kn,'^n:Art) 


exp  (j  ki  -Art)dArt  dSQ 


(22a) 


Equations  (21)  and  (22)  provide  the  equations  for  determining  the  x  and  y 
components  of  the  stochastic  transform  of  the  current  subject  to  the  neglect 
of  the  coupling  term  in  (4).  The  inclusion  of  the  coupling  term  is  not  dif¬ 
ficult  and,  in  fact,  can  be  done  by  inspection  in  view  of  its  similarity  to 
those  terms  already  included  in  (22).  Thus,  the  complete  result  is  as  fol¬ 
lows  for  q  ■  x  or  y  and  Aq  *  Ax  or  Ay; 


(2ir)  n  *> 

(2n)  n  J  J 

rx(81,k1,12.I2.-,lll.to)d6iiS2---d5n-u.  -4-  [■■■J  iy(61,l2,--,Sn 

(2w)  n  J 

rJr(61,k1.J2,J2 . SA>d*ldV-A’  2H„|c’+Jc’9k,_+)c;»k,]p,  I,, k,,-") 


1  J  *  k2x  '  k2jr|  1  1  12  2  3 

(23) 


where 


i  ff^vV 


rx(6l,kl,^2,k2’"  *  ’^n’V  “  "  (2lT)2  JJ  3Aq  [;i^k2x+:ia8^P2(kl'6o’0o'6l,k2,"^2,' 


t  ,-$  ;Ar  )exp(Jk  *Ar  )dAr  dg 
n  n  t  it  t  c 


(23a) 
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r,'B I’VVr-'-'Vv-- (2‘)2jy  Jsj  °  [j3k2/3J£2y]f2<kl-3o'V3l'k2’-V"' 

k  ;Ar  )exp(jk.  *Ar  )dAr  d8 
n  n  c  l  t  to 


(23b) 


The  z-component  of  the  stochastic  transform  of  the  current  is  obtained 
by  averaging  (3)  weighted  by  the  Fourier  kernel.  This  operation  leads  to  the 
following; 


lim 

tt+oo 


(2rr) 


1+«n 


Hi 


WV-’V-J 


3jx(Sl*V"*V 


38 


-  J 


2x 


382y  J 


*  . . .  d6lde2--,d8n  *  ° 


(24) 


Since  this  result  must  hold  for  all  values  of  ic^.k^,'-',  it  follows  that 


.Viizi 


'  3  j_— *— - +  -JL 


2x 


38 


2y 


(25) 


Discussion 

One  of  the  first  results  of  the  stochastic  transform  approach  is  illus¬ 
trated  in  (21);  that  is,  the  dependence  of  the  transformed  current  upon  the 
transverse  coordinates  x  and  y  is  identical  to  that  of  the  incident  plane  wave. 
An  important  implication  of  this  result  is  that  the  average  scattered  field 
(obtained  by  substituting  (21)  in  (8)  and  this  result  in  the  average  of  (5)) 
is  always  specular.  That  is,  the  average  scattered  field  is  non-zero  only  for 

k  -  k.  and  k  ■  k. ...  It  should  be  noted  that  the  specular  nature  of  the 
sx  ix  sy  iy 

average  scattered  field  is  a  general  and  an  exact  result  for  homogeneous  sur¬ 
face  statistics.  However,  if  the  surface  has  nonhomo geneous  statistics  then 
<Efl>  will  not,  in  general,  be  specular  because  che  joint  two-point  probability 
density  function  will  no  longer  depend  only  on  the  difference  coordinates  Ar^  . 
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The  specular  nature  of  the  average  scattered  field  was  previously  demonstrated 
by  DeSanto  [5]  using  the  diagrammatic  approach. 

It  is  interesting  to  compare  (22)*with  the  results  obtained  from  the  dia¬ 
grammatic  technique  [6].  In  the  case  of  the  latter  approach  one  obtains  an 
infinite  hierarchy  of  one -dimensional  singular  integral  equations  of  the  second 
kind.  In  order  to  obtain  the  integral  equations  in  the  hierarchy,  one  must  be 
relatively  familiar  with  the  use  of  diagram  techniques  to  form  the  appropriate 
"mass  operator"  which  is  related  to  both  the  source  term  and  the  kernel  in  the 
resulting  integral  equations.  Furthermore,  since  at  a  given  level  of  the  hier¬ 
archy,  the  resulting  integral  equaton  depends  in  part  on  the  solution  of  the 
next  lower  level  integral  equation,  the  approach  has  practical  numerical  limi¬ 
tations.  Finally,  it  is  somewhat  difficult  to  see  how  certain  statistical  fea¬ 
tures  of  the  surface  impact  the  hierarchy  of  equations.  By  contrast,  (22)  is  a 
single  integral  equation  of  the  first  kind  involving  an  infinite  number  of  vari¬ 
ables.  Because  of  the  infinite  range  of  integration,  the  integral  equation  is 
of  the  singular  type  with  respect  to  all  variables  [10,  pg.  160],  Thus,  the 
fact  that  (22)  is  an  integral  equation  of  the  first  kind  should  not  cause  too 
great  a  problem.  Obviously,  (22)  cannot  be  solved  numerically  for  much  beyond 
three  variables  (^^2)  unless  the  kernel  has  a  very  special  behavior  with 
respect  to  the  higher  dimensions.  Based  upon  these  points  of  comparison,  (22) 
exhibits  no  clear  advantages  with  respect  to  the  diagrammatic  approach  except 
that  it  is  more  compact. 

There  is,  however,  one  aspect  of  (22)  which  is  extremely  useful  and  does 
not  appear  to  be  obviously  present  in  the  diagrammatic  technique.  The  impor¬ 
tance  of  all  the  higher  order  surface  height  derivatives  is  explicitly  contained 

*The  remainder  of  this  section  will  deal  with  the  scalar  formulation  for  jq 
because  of  its  more  compact  form.  The  concepts  and  developments  apply  equally 
well  to  the  vector  equation  (23)  however. 
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in  the  two-point  joint  characteristic  function  p2(*)  in  (22a).  It  is  therefore 
a  reasonably  straightforward  task  to  determine  the  conditions  under  which  a 
certain  surface  height  derivative  is  no  longer  important  in  the  determination 
of  the  stochastic  Fourier  transform  of  the  current.  The  logic  for  accomplish¬ 
ing  this  task  is  as  follows.  Assume,  for  example,  that  one  would  like  to  know 

2 

when  the  curvature,  V  £  ,  and  all  higher  order  derivatives  have  nc^  significant 
impact  on  the  current.  If  they  are  not  important  then  J  (r  ,£,V£,  •  •  • )  will  be  essen- 

^  S  ~3  ro 

tially  independent  of  V3£, •  •  •  so  that  iq  (8j_ »^2  ’^3 ’ *  *  *  >=^q (^i* &>>  (2tt)  ^J3<5(Bn) 

because  the  transform  of  a  constant  is  a  delta  distribution.  Substituting  this 
result  in  (22)  yields; 

(2tt)3  f"'f  Jq(6l,e2)  ^1  (k1-31,k2-62,k3,---,kn)  |l 

-  r2(31,k1,32,k2,0,k3,  •  •  •  .O.k^)  j p1(k1-31,k2~B2,k3,  •  •  •  .k^)  j 

2  3 

In  order  to  make  things  explicit,  let  the  stochastic  variables  V  £,V  £,■••  be 
statistically  independent  of  X,  and  V£  .  Then  the  single-point  characteristic 
function  may  be  factored  as  follows; 

PlOtl~^l’k2~^2  ,k3,  •  ••  )  -  p1(k1-f?1,k2-32)p1(k3,k^,  •  •  * ) 
where  p^(*)  is  the  single-point  characteristic  function  for  the  indicated  vari¬ 


ables.  Substituting  this  identity  in  (26)  and  dividing  by  p^ (k3 ,k^ , • • • )  yields 


2lr)3  f'“f  Jq(Bl'62)Pl(kl"ei,k2^2)  { 


r2^rkl’®2,k2,0,k3’  ‘  *  ‘»°*kn^Aifki~ei*k2_^2)Pl(k3’kA,‘ "  ‘’kn>  }  d6lc 


2H  CH  +  j  Ch8, 

0  2  x  k- 


+  1C’SJ 


Pi(krkiZ’k2) 
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The  right  side  of  (27)  is  independent  of  whereas  the  left  side  still 

depends  on  the  variables  through  the  I^/  term.  Since  this  is  an  impossible 
situation,  one  examines  the  ^  term  to  determine  the  conditions  for  which 

it  is  essentially  (but  not  trivially)  independent  of  .  Thus .  the 

stochastic  Fourier  transform  technique  provides  a  very  orderly  procedure  for 
estimating  the  importance  of  the  surface  height  derivatives.  This  result  ap¬ 
pears  to  be  one  of  the  most  important  advantages  of  this  technique. 

Limiting  Cases 

In  the  previous  section  the  stochastic  Fourier  transform  approach  has  been 
introduced,  developed,  and  one  of  its  merits  has  been  discussed.  However,  the 
real  power  implicit  in  the  technique  results  from  its  ability  to  provide  exact 
results  for  two  very  important  limiting  cases.  In  order  to  be  specific,  the 
randomly  rough  surface  is  assumed  to  be  a  jointly  Gaussian  process;  [11]  pro¬ 
vides  details  on  the  form  of  p^  and  ^  for  a  Gaussian  process.  Such  an  assump¬ 
tion  is  not  essential  to  the  following  development  but  it  does  facilitate  com¬ 
parison  with  previous  results.  Furthermore,  it  will  be  assumed  that  the  sto¬ 
chastic  surface  height  has  a  well  defined  spatial  correlation  function 
R(Art)  =  <C(rt  )?(rt  +Art)>  which  decays  with  increasing  ]Art|  such  that  the 
concept  of  a  correlation  length  (&)  has  meaning. 

Gently  Undulating  Surface  (£>>Aq) 

The  first  limiting  case  to  be  considered  is  that  of  a  surface  of  arbitrary 
height  but  so  gently  undulating  that  none  of  the  surface  height  derivatives 
(A,  n-1,2  , •••)  are  important  to  the  determination  of  the  current.  For  such 
a  surface,  the  correlation  length  is  necessarily  much  larger  than  the  electro¬ 
magnetic  wavelength  ( l»\Q ). 

Using  the  technique  developed  in  the  previous  section  it  is  possible  to 


determine  the  conditions  under  which  Che  height  derivatives  are  not  important. 
To  shorten  matters  somewhat,  the  technique  will  be  demonstrated  by  considering 
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only  surface  slopes.  That  is,  <(Vnr;)^>  for  n*2,3,-*-  are  assumed  to  be  van  ishlngly 
small  so  that  the  single  and  two-point  joint  probability  density  functions  ap¬ 
proach  centered  delta  distributions  with  respect  to  these  variables.  Conse¬ 
quently,  the  ^oint  characteristic  functions  and  p2  in  (23)  become  indepen¬ 
dent  of  and  k  for  n»2,3,-”  and  the  transform  of  the  current  may  be  written 
n  n 

as  follows; 

s  -3  “ 

jq(61J2,?3,-*-)  =  jq(61J2)(277)  “  TT  <5(8.)  (28) 

This  result  is  merely  a  restatement  of  the  fact  that  the  current,  under  these 

2 

conditions,  is  independent  of  the  surface  curvatures  (V  £) ,  rate  of  change  of 
3 

curvatures  (V  £),  etc.  To  determine  when  the  slopes  (VQ  may  be  ignored,  is 
factored  out  of  T  ,  P  ,  and  T  in  (23),  j  (8,,30)  is  assumed  to  be  of  the  fol- 
lowing  form 

jq(3iJ2)  -  jq(B1)6(I2)(2ir)2  (29) 

and  the  residual  dependence  of  ,  r^/p^  »  and  Py/p3  uPon  ^2  examijied 

after  dividing  both  sides  of  (23)  by  p^O^).  These  operations  lead  to  the  fol¬ 
lowing  simplified  form  for  (23); 


2-n  31){1-r2(01,k1,k2)/p1(krt,k2)}d01  -  2tt  61J 

*rx(6l,k1^2)/Pl(kr01^2)dei  "  &  J  jy(31)P1(k1-B1)ry(61,k1,'k2)/p1(k1-81,k2)d31 


2H 


[cq-J<C2>k„  Cq-j<r,  2>k,  Cqlp.(k--k.  ) 
it  z  J  x  2x  x  J  y  2y  yj  1  1  iz 


(30) 


where 


^r['<l;-2>k2*+  (Bo'B1)RxJ  +  J  ifr[-  <VSy+«o-V\K} 


•  x  exp(jkt  *Art)dArtd0c 


(30a) 
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h  ’  "  fj  ^  ^  "  f<^xx]k2x+(2^-Brkl)Rx"Rxyk2x}Xe^(^it‘Att)dA^td8o 

(30b) 

^2  J f  j  ^  [<4>+Ryy]k2y+(2V^-kl)Ry-Rxyk2y}Xexp(jkit*A?t)dA^td61 


Pi  (2TV)' 
and 

x  .  exp  j-  [<?2>-RXB0-81XB0-k1)  +  (S^Xk^R,  +k2yRy)} 


o 

(30c) 

(30d) 


2  2 

In  the  above  expressions  <Z  >  and  >  are  the  variances  of  the  slopes  in  the 

x  y 

x  and  y  directions  while  Rx=3R/3Ax  ,  Ry*3R/3Ay  ,  etc. 

Consistent  with  the  assumption  in  (29),  j^(B^)  is  to  be  determined  from 
(30)  and  the  resulting  solution  must  be  independent  of  the  slope  statistics 

•f 

and  k^.  For  arbitrary  slope  variances  this  is  not  possible  and  it  is  necessary 
to  examine  the  conditions  under  which  the  slope  dependent  terms  in  (30)  can  be 

made  small.  The  terms  on  the  right  hand  side  of  (30)  which  depend  linearly 

-*•  2  2 

upon  k„  can  be  made  small  by  taking  <£  >  <<  1  and  <£  >  <<  1  .  If  the  slope 

x  y 

variances  are  small,  the  correlation  function  may  be  approximated  as  follows; 


->  2  12  212  2 
R(Art)  =  <C  >  -  |  <Sx>  (Ax/  -  y<Cy  >(Ay) 


(31) 


,n_N2, 


because  < (V  £)  >  for  n=2,3,'**  have  already  been  assumed  to  be  vanishingly 
small.  Also,  the  coordinate  system  has  been  oriented  such  that  there  are  no 
Ax* Ay  terms  in  (31).  In  view  of  (31),  the  terms  in  (30b)  and  (30c)  which  de¬ 
pend  linearly  on  k2x  and  k2y,  respectively,  vanish.  The  exponential  in  (30d), 
using  (31),  may  be  written  as  follows; 

X  -  exp  |  -  |[<CX2>  (Ax)2  +  <?y2>  (Ay)2](B0-B1)(Bo-k1)| 


2 


n«o 


(-Dn  (Wn  [<  ^x2>k2xAx  +  <  S2>k2VAy1 n 


n! 


(32) 


The  terms  in  (30a)- (30c)  which  have  R  or  R  multiplying  the  Green's  function 

x  y 
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derivatives  may  be  lumped  with  the  series  expansion  in  (32).  The  Ar  -integration 

->■ 

in  (30a)- (30c)  may  now  be  treated  as  a  convolution  of  the  Ar^-Fourier  transform 
of  (32)  with  the  transform  of  the  Green's  function  derivatives.  The  transform 
of  the  Gaussian  factor  in  (32)  is 


2Tr(eo-3i)(eo-ki)/<Cx2><Cy2> 


2<W(3o-ki)  L 


? 

k  y 

<s2> 


For  small  slope  variances  (32a)  behaves  as  a  delta  distribution  which  combs  the 
transform  of  the  Green's  function  derivatives  out  of  convolution.  The  terms 
in  the  power  series  in  (32)  give  rise  to  n^  derivatives  of  the  delta  distribu¬ 
tion  which,  in  turn,  yields  derivatives  of  the  transform  of  the  Green's  func¬ 
tions.  However,  because  of  the  slope  variances  in  the  power  series  in  (32), 
the  n^h  derivative  of  the  transform  of  the  Green's  functions  will  be  multiplied 
by  a  slope  variance  raised  to  the  nth  power  which  is  vanishingly  small.  Thus, 
the  only  term  left  in  (30a)- (30d)  is 


■>■*// 


8^  exp(jkit*Art)  dArtdpo 


However  since  the  Ar  -Fourier  transform  of  £  is  an  odd  function  of  6  ,  T,/p 

t  £  O  c  1 

is  identically  zero  and  (30)  becomes 


2?/  yei>5i<krBi,dsi  =  y’yw  »3> 

so  that  j{J(81)  =  2HoCzq  6(6L-klz)(2TT)  or 

jq(z)  =  2Hq  C  q  exp(-jkizC)  (34) 

which  is  the  scalar  single  scattering  physical  optics  result.  With  = kgz  = -kiz , 

(33)  shows  that  the  average  field  scattered  by  the  surface  is  attenuated  from 

~  2  2  2 

its  flat  plane  value  by  the  factor  Pi(-2kiz)  *  exp(-2kQ  >cos  0).  This  is 
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the  classical  result  which  simply  says  that  if  the  surface  is  sufficiently 
smooth,  the  average  scattered  field  appears  to  be  due  to  a  randomly  elevated 
plane.  1  he  important  point  to  note  from  this  development  is  that  for  a  very 
gently  undulating  surface  which  has  a  correlation  length  which  is  large  com¬ 
pared  to  Xq,  scalar  single  scattering  or  physical  optics  is  applicable. 

Uniformly  Rough  Surface  (£  <<  X^)^ 

A  much  more  interesting  limiting  case  is  provided  by  surfaces  which  have 
arbitrary  height  but  a  correlation  length  which  is  small  compared  to  X  .  Such 
surfaces  necessarily  have  a  great  deal  of  their  height  variance  contributed  by 
undulation  frequencies  the  order  of  and  greater  than  k^ ,  the  electromagnetic 
wavenumber.  Consequently,  there  is  absolutely  no  reason  to  expect  physical 
optics  to  be  an  adequate  approximation  for  the  scattering  process  and,  indeed, 
this  is  found  to  be  the  case. 

Before  proceeding,  it  is  worthwhile  pointing  out  why  it  is  the  ratio  of 

the  correlation  length  to  Xq  that  is  a  critical  parameter.  Departures  from 

single  scattering  are  caused  by  the  ,  T  ,  and  T  terms  in  (23).  These  terms, 

/  x  y 

as  illustrated  in  (22a),  (23a),  and  (23b),  are  determined  by  an  integration  over 
-► 

Arfc.  The  quantities  which  appear  in  the  Art~integrands  are  dependent  upon  Xq 
and  the  correlation  length  as  scale  parameters.  Thus,  these  parameters  deter¬ 
mine  whether  it  is  the  electromagnetic  functions  or  the  statistical  functions 
which  dominate  the  Art~ integration. 

When  the  correlation  length  is  much  smaller  than  X  ,  the  statistical  func- 

o 

tions  in  (23)  vary  much  more  rapidly  with  Arfc  than  the  electromagnetic  functions. 
In  fact,  the  statistical  functions  go  from  completely  correlated  to  completely 


•it  is  difficult  to  describe  a  surface  having  l  «  Xo  with  a  few  short  words. 
"Uniformly  rough"  is  appropriate  when  the  correlation  function  is  a  Gaussian 
function  of  the  form  exp (-] | ^/t?) ,  in  this  case  the  surface  height  spectrum 
is  also  Gaussian  and  if  &  «  Xr 


then  the  roughness  is  spread  nearly  uniformly 


throughout  all  spatial  frequencies  kQ 


2tt/X, 


o’ 
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decorrelated  over  a  region  of  essentially  zero  measure  when  compared  to  the 
scale  of  variation  of  the  electromagnetic  functions.  Consequently,  the  statis¬ 
tical  functions  in  (23)  may  be  replaced  by  their  values  for  complete  decorre¬ 
lation,  i.e.  |Art|  -*•  00 ,  with  no  appreciable  error.  This  implies  that  the  cor¬ 
relation  function  and  all  of  its  higher  order  derivatives  are  set  to  zero  in 
(22a),  (23a),  and (23b).  Using  this  approximation  and  the  fact  that  decorrela¬ 
tion  implies  statistical  independence  fcr  a  Gaussian  process,  p£  in  the  expres¬ 
sions  for  r«,  T  ,  and  T  simplifies  to  the  following  form  [11]; 
z  x  y 


=  exp  [- <-2>(Bo-6i)  (6o_ki)]  Px (k1~ 3 x ) 


•  JT2  pA>Pi(A> 


where  P^(k^)  and  p^(-i^)  are  the  characteristic  functions  for  V1*;.  Substituting 


these  simplifications  in  (23)  and  noting  that 


Pl(ki-gi,k2-|2.  "  ‘  .kQ-ln)  “  PiOtj^i)  i^O^-j^)  » 

i=2 

it  is  found  that  (3^,^ > » Bn)  has  the  following  dependence  upon  ^i>i=3,4,-*-  ; 

s  -3  n 

j(]L  1)  =  j  (3.  ,P9)  (2ir)  n  TT  i(B.)  (35) 

q  1  z  n  q  i  L  x 

n  _  -+■ 

Using  this  in  (23)  and  dividing  out  the  common  dependence  upon  'TT  p.(k.)  results 

i=2  1  1 

in  the  following  two  coupled  equations  for  the  current  components; 

j/vVVPi0 -V  {v^l’V  exp  [^x2>k2xB2x  +  >k2y^2y  ] 


1’^  J  D 

[J 


<?xHx^x-j  <S2>k2yS+ljxd3o 


<Cx  >(k2y  +  e2y)  ®x 


]  X  d30 J  dexc 


2H  fc  x  -  j<£  2>k„  Cx]p.(k1-k,  ) 
oL  x  J  y  2y  y  J  *1  J  iz 


(36a) 
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JJ  (vw  **>•  K2 

-yBA>  J  [-j< 


>k2xB2x  *  VV2y] 


2>k9„i,  +  j  <ry>e‘2y*y  +  sJxdPD 


-jx(P 


I'V J  [^^><k2x+e2x>S]  *  d6o}  dBldB2  *  2Ho[CzV-J^x>k2xCx]i’l<krklz) 


(16b) 


where  the  triple  tilde  symbol  over  the  Green's  function  derivatives  denotes  the 
three-dimensional  Fourier  transform  with  respect  to  Ax,  Ay,  and  A r„  The  func¬ 
tion  x  given  by 

X  -  (2rr)~2  exp  [ -<?2>(go-ei)(g0-k1)]  . 

Expanding  the  exponential  factor  containing  k2x  and  k^  in  a  series  and  equating 
like  powers  of  k2x  and  k^  on  both  sides  of  (36),  it  becomes  clear  that 
jq(61J2)  may  be  expressed  as  follows; 

WV-  {jx1>«V4«2*>6<82y>  +  3x^  (By3B(B2x34'  y  ^  ^  ^ 


dy(8l’^2>=  <81>4 (B2x^B2yk  +  3y2> <S1>«'  «2x,6<62y)} <2"’ 


(37) 


Substituting  these  results  in  (36)  and  equating  like  coefficients  of  k9^  and 
k^y  yields  the  following  four  equations  for  the  8^-dependent  parts  of  the  trans¬ 
formed  current; 

s  J  i(w{‘-^5  / ir."i>[-<^(V8i><Vki,]<l6o)',er2Hoc°Pi<krklz) 


(38) 


5  /  Pi<krt>0r<V^^,(V  Jv  -xp[-<c2><eo-0iM6o-k13j  de0 


(2tt) 


_  An 

2  V 


•>tsi) J  *,'^[-<^><VBi>(Bo-ki>]dBo}d0i  2Hoc,q’  h(krkix)  <391 
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where  q  and  q'  are  paired  as  q*x,  q'  =y  and  q  =  y,q'  =  x.  Using  (37)  in  (21), 
(5),  and  (8)  shows  that  the  part  of  the  current  which  depends  upon  6  *  C3 2 q ^  con_ 
tributes  nothing  to  the  average  scattered  field  which  from  (5)  may  be  expressed 
as  follows; 

<E  (k  )>  =  jirk  n  g(R  )<5(k  -k.  )[k  xf  (k  )]  (4( 

s  sz  000  sis  s  q  sz 


where  for  q*x  or  y. 


T  (k  )  =  ~ 
q  sz  2tt 


J  ■iq1><kl>5l  <ksz-kl)dkl 


and  ^sz=~^£z  because  of  the  specularity  of  the  average  scattered  field.  Sub¬ 
stituting  (41)  in  (38),  combining  p^(k^-B^)  with  the  exponential  term  [11],  and 


noting  that 


g  (£  )  =  j  4 it  lim  -j - | - 

5  e-K)  8  -  k.  +  j 


yields  the  following  integral  equation  for  T^(k)  ; 


T  (k)  -  2H 
q  o 


[q*  (z  x h)  ]p.  (k-k  )  +  i  J  limf-r - ^ - -1  p. (k-B.)T  (Bn)dB 

1  iz  *  J  e-+o|_B  -k.2  +jej  0  q  0  o 


Both.  (42)  and  some  of  the  interim  results  leading  to  it  have  significant 
implications.  Eqns.  (35)  and  (37)  imply  that  the  product  of  the  surface  cur- 

/  2  T 

rent  and  /I  + 1,  +  C,  may  be  divided  into  two  parts  with  one  part  totally  inde- 

x  y 

pendent  of  the  surface  height  derivatives,  Vn£  .  It  is  tempting  to  declare 

that  the  other  part  depends  linearly  on  one  of  the  slope  components  since  this 

would  appear  to  follow  from  the  o'  (8„  )  in  (37)  and  the  Fourier  transform  rcla- 

?.q 

tionship  between  and  j  ,  see  (11)  and  (21).  However,  it  is  not  clear  that 


such  a  Fourier  inversion  unique 


recovers  the  true  stochastic  nature  of  the 


current,  i.e.  the  true  dependence  on  rt  .  This  simply  says  that  one  cannot, 
in  general,  go  from  j  back  to  Jg*q  • 

Using  (38)  and  (39)  and  the  fact  that  C  X =  -  C ^  ,  it  can  be  shown  that 

y  ^ 

J*Ji2>(Bi)?i(krBi)dBi  =“ /jJ2)(Wkr*i)dBi 


2-26 


(2) 

and  that  j  ((3^)  can  be  determined  from  knowledge  of  Tx(k)  according  to  the 
following; 


H  J  Ji2>cei>Pi(ki-ei)d0i 


■  -  j^»ocy  h(krki2>  +|ik...-k 


ty 


/ 


IX  cx 

z 


lim  /  — y - 2 - 

e+0  )  3  “  -  k  f  +  jc 

'  o  iz 


P.  (k,-f?  )T  (B  )d6 

1  1  O  X  o  o 


(44) 


x  (2 ) 

for  C  ^  0,  It  is  important  to  realize  however  that  j  contributes  nothing 

Z  q  (2) 

to  the  average  scattered  field.  Discussion  of  the  effect  of  i  on  the  mean 

q 

square  scattered  field  will  be  left  to  future  investigations. 

Equation  (42)  is  very  similar  to  the  integral  equation  obtained  by  DeSanto 

[6]  for  this  "uniformly  rough"  surface  limit.  It  differs  in  that  DeSanto  ob¬ 
tains  a  in  the  denominator  of  the  kernel  of  the  integral  term  rather  than 


in  the  numerator  as  in  (42).  However,  DeSanto's  integral  equation  results  from 


only  the  first  term  in  his  diagrammatic  hierarchy  of  equations  whereas  (42)  is 

exact.  It  is  therefore  not  surprising  that  the  two  integral  equations  differ. 

At  the  beginning  of  this  section,  it  was  stated  that  the  gently  undulating 

and  uniformly  rough  surfaces  are  very  important  limiting  cases.  Although  it  is 

difficult  to  prove  mathematically,  there  are  strong  physical  agruments  to  suggest 

that  these  two  cases  provide  lower  and  upper  bounds  on  the  average  scattered 

2 

field.  That  is,  for  a  fixed  height  variance,  >,  the  gently  undulating  sur¬ 


face  gives  a  lower  bound  on  the  coherent  field  while  the  uniformly  rough  surface 
yields  the  upper  bound.  As  noted  previously,  the  gently  undulating  surface  be¬ 
haves  as  a  randomly  eleva’ed  plane  which,  in  turn,  puts  all  of  the  randomness 
of  the  surface  in  the  phase  of  the  scattered  field.  Upon  averaging  over  an  en¬ 
semble  of  scattered  fields,  the  randomness  of  the  phase  leads  to  a  maximum 

effect  on  the  mean  scattered  field  and,  consequently,  a  minimum  value.  As  the 
surface  approaches  the  uniformly  rough  case,  the  randomness  of  the  surface  no 
longer  appears  totally  in  the  phase  of  the  scattered  field  due  to  multiple 
scattering  on  the  surface.  This  is  just  a  restatement  of  the  fact  that  the  phy¬ 
sical  optics  approximation  is  not  accurate  for  surface  features  having  an  undu¬ 


lation  frequency  which  is  small  compared  to  k  .  Thus,  once  the  correlation 
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length  of  the  surface  exceeds  the  electromagnetic  wavelength,  the  effects  of 
the  random  height  on  the  phase  of  the  scattered  field  are  minimized.  Conse¬ 
quently,  this  should  lead  to  a  maximum  coherent  scattered  field. 

One  final  interesting  facet  of  these  two  limiting  cases  is  the  points  which 
they  share  in  common.  Both  cases  depend  only  on  the  3g/<HQ  term  in  the  integral 
in  (10).  For  the  gently  undulating  surface  the  contribution  from  this  term  is 
essentially  zero.  For  the  uniformly  rough  surface  this  term  encompasses  all 
of  the  multiple  scattering  that  is  important  to  the  coherent  or  average  scattered 
field.  A  second  point  of  commonality  is  the  fact  that  the  currents  which  con¬ 
tribute  to  the  average  scattered  fieldfare  independent  of  all  the  surface  height 
derivatives.  These  facts  form  the  basis  for  an  approximate  approach  to  solving 
(42)  which  is  detailed  elsewhere  [12]. 

®For  the  gently  undulating  surface,  this  is  the  total  current  while  for  the 

uniformly  rough  surface  it  is  only  a  portion  of  the  current  which  contributes 

to  <t  >• 
s 
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Abstract 

An  approximate  solution  for  the  average  field  scattered  by  a  perfectly 
conducting  randomly  rough  surface  having  a  correlation  length  much  smaller  than 
the  electromagnetic  wavelength  is  presented.  The  analysis  is  based  on  the  use 
of  a  substitute  surface  which  gives  rise  to  the  same  describing  equations  as 
the  true  surface  relative  to  the  average  scattered  field.  The  substitute  sur¬ 
face  comprises  large,  non-overlapping,  flat  areas  having  random  elevations 
with  respect  to  the  mean  planar  surface  and  arbitrary  correlation  between 
adjacent  areas.  The  average  scattered  field  is  shown  to  depend  upon  the  number 
of  interacting  areas  and  the  surface  roughness.  For  a  given  range  of  surface 
roughness  there  is  a  specific  number  of  interacting  areas  which  dominate  the 
average  scattered  field.  It  is  demonstrated  how  this  number  can  be  computed 
and  how  a  continuous  curve  of  average  scattered  field  as  a  function  of  surface 
roughness  is  obtained.  Of  particular  importance  is  the  quantitative  correspon¬ 
dence  established  in  this  paper  between  the  surface  roughness  and  the  degree 
of  multiple  interaction  on  the  rough  surface. 

Introduction  and  Summary 

Single  scattering  approximations  have  served  very  nicely  to  illustrate 
some  of  the  salient  features  of  scattering  from  randomly  rough  surfaces  [1]. 
There  is,  however,  a  lingering  controversy  as  to  how  far  these  approximations 
can  be  extended  [2].  The  basic  problem  associated  with  establishing  the  range 
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of  validity  of  these  approximations  is  that  the  next  order  correction  is  very 
difficult  to  obtain.  A  number  of  solutions  have  been  presented  which  are  for¬ 
mally  exact  but  are  practically  limited  to  a  two-term  perturbation  approxima¬ 
tion  [3,4].  The  pioneering  work  of  Zipfel  and  DeSanto  [5]  provided  the  first 
approximate  correction  to  single  scattering  theory  which  did  not  require  a 
small  perturbation  in  the  surface  height.  Subsequent  numerical  work  by  DeSanto 
and  Shisha  [6]  showed  that,  for  a  Gaussian  distributed  surface,  single  scat¬ 
tering  grossly  underestimated  the  strength  of  the  average  scattered  field  for 
large  surface  roughness  and  small  correlation  length.  Unfortunately,  the 
detailed  mathematics  essential  to  obtaining  this  correction  to  the  single 
scattering  approximation  did  not  lend  itself  to  a  straightforward  physical 
interpretation  [7]. 

Recently  another  approach  to  the  rough  surface  scattering  problem  has 
been  developed  [8].  Apart  from  starting  with  the  magnetic  field  integral 
equation,  the  unique  feature  of  this  approach  was  the  conversion  of  averages 
to  convolutions  in  stochastic  transform  space.  The  purpose  of  this  new  approach 
was  not  to  duplicate  the  excellent  work  of  Zipfel  and  DeSanto  but  rather  to 
obtain  a  solution  which  clearly  illustrated  the  range  of  validity  of  approxi¬ 
mations  necessary  to  simplify  the  full-blown  problem.  One  of  the  results  of 
this  effort  was  the  derivation  of  an  exact  compact  integral  equation  for  the 
average  scattered  field  in  the  case  of  a  surface  having  a  correlation  length 
much  smaller  than  the  electromagnetic  wavelength.  This  equation  was  similar 
to  the  approximate  integral  equation  obtained  by  Zipfel  and  DeSanto  [5].  In 
[8]  it  was  shown  that  the  average  field  scattered  from  a  randomly  rough  per¬ 
fectly  conducting  surface  was  attenuated  from  its  flat  surface  value  by  the 
quantity  T(k=*k  *>  -  k  )  where  k  is  the  z-component  of  the  vector  wave- 
number  in  the  scatter  direction  and  k  is  the  z-component  of  the  vector 
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wavenumber  in  the  incidence  direction.  The  condition  k  «- k,  results  from  the 

sz  i  z 

specularity  of  the  average  scattered  field.  For  Gaussian  surfaces  having  a  corre¬ 


lation  length  (£.)  much  smaller  than  the  electromagnetic  wavelength  (Aq),  the 
exact  equation  for  T^(k)  was  found  to  be  the  following; 


Tq(k)  =  2HQ[q-(z  *h)]  PL  (k-kiz)  +  i 


+  ?  f  nmT-^A - 

J  c-*°k  -ki2+^ 


p  (k-6  )T  (B  )dB  (1) 
1  o  q  o  o 


where  and  h  are,  respectively,  the  amplitude  and  polarization  of  the  mag¬ 
netic  field  incident  upon  the  surface.  In  (1)  q  is  either  x  or  y  and  is  the 
Fourier  transform  of  the  surface  height  probability  density  function. 

A  slightly  different  form  of  this  equation  (with  3  in  the  denominator 

o 

of  the  bracketed  term  rather  than  in  the  numerator)  was  solved  numerically  by 
DeSanto  and  Shisha  [6]  and  it  would  appear  to  be  a  straightforward  task  to 
apply  their  technique  to  (1).  However,  this  writer  has  some  concerns  about 
using  the  same  technique  to  numerically  solve  (1).  First,  although  the  dif¬ 
ference  between  (1)  and  the  equation  solved  in  [6]  may  appear  to  be  minor,  it 
is  likely  that  the  difference  may  significantly  alter  the  convergence  properties 
of  the  integral  in  (1).  This,  in  turn,  leads  to  the  distinct  possibility  of 
singularities  in  T  (k)  for  k  t  -  k  which  could  not  be  accounted  for  using  the 

4  1 Z 

quadratic  spline  technique  employed  in  [6].  In  short,  the  numerical  solution 
of  (1)  may  require  a  great  deal  more  in-depth  analysis  than  presented  in  [6]. 

Finally,  there  is  a  very  practical  limitation  contained  in  a  numerical 
solution  of  (1).  Assuming  that  a  suitable  technique  can  be  devised  for  solving 
(1)  numerically,  what  one  ends  up  with  are  numbers!  While  these  numbers  are 
certainly  important,  they  add  nothing  to  our  understanding  of  the  complex 
process  of  multiple  scattering.  If  real  progress  is  to  be  made  in  moving 
beyond  our  present  knowledge,  it  must  come  from  a  clear  understanding  of  the 
physics  of  the  problem  and  not  just  numbers.  Consequently,  considering  the 
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potential  difficulties  associated  with  a  numerical  solution  of  (1)  along  with 
the  limited  insight  available  from  such  results,  there  is  a  definite  need  for 
an  alternate  solution  of  this  problem.  The  purpose  of  this  paper  is  to  present 
one  such  alternate  approximate  solution. 

The  cornerstone  of  the  new  solution  is  the  previously  derived  [8]  similar¬ 
ity  between  the  stochastic  Fourier  transform  of  the  total  current  induced  on 
a  gently  undulating  surface  (4  »  X  )  and  the  part  of  the  stochastic  Fourier 
transform  of  the  current  induced  on  a  uniformly  rough  surface  ( £  <<  X  )  which 
is  responsible  for  the  coherent  scattered  field.  That  is,  both  of  the  trans¬ 
formed  quantities  have  the  same  dependence  on  all  surface  height  derivative 
statistics  and  they  satisfy  exactly  the  same  reduced  form  of  the  transformed 
magnetic  field  integral  equation.  This  similarity  suggests  that  if  the  gently 
undulating  surface  ( £.  »  Xq)  could  be  altered  so  that,  from  a  mathematical  point 
of  view,  the  case  of  SL  «  XQ  could  be  addressed  then  this  altered  surface  could 
be  substituted  for  the  uniformly  rough  surface  .  Furthermore,  if  this  substi¬ 
tute  surface  is  such  that  its  scattering  properties  with  respect  to  the  aver¬ 
age  field  could  be  determined  then  these  results  would  apply  also  to  the  uni¬ 
formly  rough  surface.  The  key  to  success  in  such  an  approach  is  assuring  that 
the  substitute  problem  can  be  solved  more  easily  than  the  original  problem 
and  that  the  results  are  in  complete  agreement  with  the  physics  of  the  problem. 

The  transformation  of  the  gently  undulating  surface  (Z  >>  X  )  into  a  sub¬ 
stitute  surface  for  the  uniformly  rough  surface  (2.  <<  X  )  proceeds  as  follows. 
Because  of  the  smallness  of  the  surface  slopes,  curvatures,  etc.  ,  the  gently 
undulating  surface  may  be  approximated  by  an  infinite  number  of  large 

*It  is  assumed  that  the  original  similarity  between  the  transformed  currents 
is  invariant  under  the  alteration  process. 
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non-overlapping  conducting  areas  which  are  flat  (non-inclined)  but  randomly 
elevated  with  respect  to  the  z*0  plane.  Toe  current  induced  on  any  one  flat 
area  is,  however,  determined  by  only  a  finite  (but  unknown)  number  of  neigh¬ 
boring  flat  areas.  An  exact  determination  of  the  current  on  the  rough  surface, 
in  this  approximation,  is  possible  and  it  is  found  to  depend  upon  the  relative 
elevations  of  the  non-overlapping  areas  in  a  rather  simple  and  physically 
meaningful  manner.  The  complex  amplitude  of  the  average  scattered  field  is 
obtained  by  averaging  the  product  of  the  random  height  dependent  part  of  the 

current  and  the  Fourier  kernel  exp(j  k  £)  where  k  is  the  z-component  of  the 

sz  s  z 

wavenumber  vector  in  the  scattering  direction  and  C,  is  the  random  surface  height. 

In  order  to  accomplish  the  field  averaging  one  must  have  knowledge  of  the 
joint  probability  density  function  for  the  heights  of  the  electromagnetically 
important  neighboring  areas.  This,  in  turn,  raises  the  question  of  the  degree 
of  correlation  between  these  areas.  Thus,  it  is  at  this  point  that  total 
decorrelation  may  be  inserted  to  complete  the  similitude  between  the  substitute 
surface  and  uniformly  rough  surface,  i.e.  since  Z  «  \q  the  uniformly  rough  sur¬ 
face  may  be  approximately  characterized  by  Z  ~  0  or  complete  decorrelation. 

The  substitute  surface  therefore  comprises  an  infinite  number  of  large  non¬ 
overlapping  conducting  areas  which  are  flat  but  randomly  elevated  and  uncor¬ 
related  with  each  other. 

In  the  case  of  surfaces  for  which  decorrelation  implies  statistical  inde¬ 
pendence  (such  as  a  Gaussian)  the  final  result  reduces  to  one  integration  which 
can  be  accomplished  numerically.  Not  surprisingly,  the  average  scattered  field 
is  found  to  depend  upon  the  number  of  areas  or  regions  on  the  rough  surface 
allowed  to  interact  with  the  point  in  question.  For  a  fixed  number  of  inter¬ 
acting  areas,  there  is  a  specific  range  of  roughness  or  Rayleigh  parameter 
for  which  the  average  scattered  field  is  a  maximum.  As  the  roughness  decreases 
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toward  zero,  the  number  of  interacting  areas  on  the  surface  which  maximize 
the  scattered  field  also  approaches  zero.  Conversely,  as  the  roughness  in¬ 
creases  the  number  of  interacting  areas  necessary  to  maximize  the  average 
scattered  field  increases.  In  terms  of  the  physics  of  the  scattering  process, 
these  results  have  a  very  clear  meaning.  In  essence  the  solution  comprises 
a  discrete  approximation  to  the  fact  that,  as  the  surface  roughness  increases, 
the  current  at  any  point  on  the  surface  is  dependent  upon  an  increasing  area 
of  the  surrounding  surface  due  to  multiple  interaction  and  not  just  the  behavior 
of  the  surface  in  the  immediate  vicinity  of  the  point  in  question.  This  result 
represents  progress  in  the  understanding  of  scattering  from  rough  sur¬ 

faces  in  that  it  is  the  first  time  that  the  complicated  mathematics  of  mutual 
interaction  have  been  simplified  to  the  point  where  the  result  can  be  put  into 
one-to-one  correspondence  with  the  physics  of  the  process. 

Numerical  results  for  the  strength  of  the  average  or  coherent  scattered 
field  are  presented  for  the  case  of  a  Gaussian  distributed  surface.  These 
results  show  a  significantly  stronger  coherent  power  for  large  surface  rough¬ 
ness  or  Rayleigh  parameter  than  is  predicted  by  single  scatter  or  the  physical 
optics  approximation.  The  apparent  reason  for  the  failure  of  the  single  scat¬ 
ter  approximation  is  as  follows.  The  single  scatter  or  physical  optics  approxi¬ 
mation  attributes  the  effect  of  the  surface  roughness  to  a  phase  modification 
which,  when  averaged,  tends  to  maximize  the  effect  of  the  roughness.  For 

SL  «  \  ,  the  solution  obtained  in  this  paper  shows  that  mutual  interaction 
o 

gives  rise  to  the  same  phase  modification  as  obtained  with  single  scatter  plus 
an  additional  amplitude  effect.  However,  this  amplitude  effect  deemphasizes 
the  influence  of  the  surface  roughness  as  predicted  by  the  phase  only  approxi¬ 
mation.  This  result  also  provides  some  insight  into  why  iterative  or  "back  and 
forth"  multiple  scattering  attempts  to  solve  the  problem  are  successful  only 
if  the  iterative  series  can  at  least  be  partially  summed.  That  is,  such  attempts 
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lead  to  series  terms  which  contain,  primarily,  a  phase  perturbation;  thus,  no 


single  term  or  any  finite  sum  of  terms  can  possibly  describe  the  important 
amplitudes  effects  found  in  this  paper. 

A  comparison  with  the  results  obtained  by  DeSanto  and  Shisha  [6]  is  also 
presented.  While  there  is  agreement  in  predicting  a  greater  average  scattered 
field  than  single  scattering  for  large  surface  roughness,  there  is  also  a  dif¬ 
ference  in  the  detailed  behavior.  This  is  not  surprising  given  that  the 
DeSanto  and  Shisha  result  stems  from  only  the  first  integral  equation  in  an 
infinite  hierarchy  of  integral  equations. 

Determination  of  the  Surface  Current 

The  rough  surface  provides  the  boundary  between  free-space  and  a  perfectly 
conducting  medium  and  it  is  assumed  to  be  infinite  in  extent  in  the  x  and  y 
coordinates  of  a  conventional  (x,y,z)  coordinate  system.  The  random  surface 
height  C(x,y)  measured  from  the  z  *»  0  plane  is  assumed  to  comprise  a  zero 
mean  statistically  homogeneous  process  with  the  mean  surface  equal  to  the  z  =  0 
plane.  The  magnetic  field  integral  equation  110,  pg.  354]  for  the  current 

-►  a  a  /s  .  ... 

J  (r)  induced  on  the  surface  at  a  point  r  =  xx  +  yy  +  £z  due  to  an  incident 

8 

"*i  -+  /s  -► 

magnetic  field  of  the  form  H  (r)  =  HQh  exp  (- j  k^*  r )  may  be  expressed  in  the 
following  form  [8]; 

J(r)  *  2N(r)*Hi(r)  +  j  [N(r)  *Vog)  J  (rQ)  -  [N(r)*  J  (rQ)  ]  VQg  |  drt  (2) 

where  J(r)  is  related  to  Jg(r)  by 


J(£) 


/ 


1  +  Oc/Sxr  +  OC/9y)' 


J  (r) 
s  ' 
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and 


>  +  3  r  3  cj  A 

N(r)  *  -  ^Jx  “g^y  +  z  .  The  Incident  wavenumber  vector 
kixx  +  klyy  +  k±z  z  is  given  by 


k^  =  -  kQ  [sin0  cos4>  x  +  sinB  sin0  y  +  cos0  z] 


where  0  is  the  angle  of  incidence  with  respect  to  the  normal  to  the  mean  sur¬ 
face,  <J>  is  the  azimuth  angle  measured  from  the  x-axis,  and  k  =  2tt/A  is 

o  o 

the  free-space  wavenumber.  The  quantities  3£/3x  and  3c;/3y  are 

the  x  and  y  components  of  the  surface  slopes  at  the  point  r  on  the  surface, 

r  =  x^x  +  yQy  ,  and  g(jr  -  r  |)  is  proportional  to  the  free-space  Green's 
o 

function,  i.e. 


exp(  -  j  kjr-  rj) 


(3) 


Using  the  fact  that  the  current  must  be  orthogonal  to  the  surface  normal,  i.e. 
-►  ■+ 

N(ro)  •  J( r^)  «  0  ,  leads  to  the  following  relationship  between  the  vector  com¬ 
ponents  of  the  current; 


J  (r  )  -  (3C  /3x  )  J  (r  )  +  Or  /3y  )  J  (r  ) 

Z  O  O  O  X  O  O  o  y  o 


(4) 


Using  (4),  it  can  be  shown  that  the  term  [N(r)  *J(ro>]  in  (2)  is  proportional 
to  the  difference  between  the  surface  slopes  at  the  points  r  and  rQ  on 
the  surface.  This  term  contributes  nothing  to  that  part  of  the  current  respon¬ 
sible  for  the  coherent  scattered  field  for  either  £  »  A  or  t  <<  A  [8}  and 

o  o 

so  (2)  reduces  to  three  scalar  equations  or,  alternately,  (4)  and  the  following 
two  expressions; 


\(h 


2^- [N(r)  x  Hi(r) ]  +  ^ 


3 K  _3g_  3C  3g  , 

3x  3x  "  5y  3y  + 
o  o 


J 

q 


(r  )dr 
o'  t 

o 


(5) 
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where  q  ■  x  or  y  and  q  is  a  unit  vector  in  either  of  these  directions. 
The  slope  dependent  terms  in  (5)  also  can  be  shown  to  contribute  nothing  for 
either  case  [8]  and  (5)  becomes 


Jq(r) 


JqUr) 


2n 


f 


_l£_ 

J  no 


J  (r  )  dr 
q  o  t 


(6) 


for  q  *  x  or  y  and  J^Cr)  =  2q*(z  *  H^Cr)} 

The  argument  of  the  Green's  function  in  (6)  is  ✓iV;,  +  <^-V2  ■ 

o 

Starting  with  the  gently  undulating  surface  (£>>Aq)  and  using  the  fact  that 

the  surface  curvatures  and  slopes  are  very  small  implies  that  the  quantity 

(£  -  £  )  is  nearly  a  constant  over  a  large  range  of  r..  .  Thus,  (6)  may  be 

o  lo 

approximated  as  follows; 


J  (r)  ~  J 
q 


i 

q 


3(C  -Cn) 


J  (r  )  dr 
q  n 


t 

n 


(7) 


n 

where  A  is  the  area  of  the  z  *  0  plane  over  which  £  -  £  is  nearly  con- 
n  n 

stant  and  M+l  is  the  number  of  these  areas  which  make  up  the  entire  z  =  0 
plane.  To  include  the  entire  z  ■  0  plane  in  the  integration  in  (7),  M  should 
be  infinite;  however  this  ignores  the  fact  that  only  a  finite  area  about  the 
point  (rt,C)  is  effective  in  determining  Jq(rt,£)  *  Consequently,  M  will  be 
replaced  by  N  which  is  assumed  to  be  finite. 

In  essence  (7)  is  a  discrete  approximation  to  (6);  retention  of  the  inte¬ 
gration  in  (7)  merely  reflects  the  fact  that  (C  -  C  )  is  constant  to  some 

prescribed  tolerance  over  an  area  An  and  not  just  a  point.  Defining  sn(rt  ) 

11  n 

to  be  the  support  function  for  the  nth  area,  i.e.  S  =1  for  r  £  A  and 

n  cn  n 

zero  otherwise,  (7)  may  be  rewritten  as  follows; 


Jq(r) 


OC 

n  r 

J 

n«l  -oo 


38(rt-r  .C-?n> 

_ n _ 


W  Sn  (rt  )drt 


(8) 
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r 


where  C-C  is  assumed  to  be  constant  or  independent  of  r  .  Figure  1 
n  ln 

illustrates  what  has  been  done  to  the  surface  in  going  from  (6)  to  (7).  In 

essence  the  continuous  surface  has  been  approximated  by  a  non-overlapping  col¬ 
lection  of  discrete  flat  areas  which  are  all  parallel  to  the  z  =  0  plane. 

The  derivative  of  the  Green's  function  may  be  expressed  in  terms  of  a 
two  dimensional  inverse  Fourier  transform  {5]; 


t  ,Kn)  r  _ i 

A? -  "  2  J  sgn  (A^a)eXp  ['■ j  ^ ko”  I  ^  j 


3g(Art  ,AcJ 
_ n 

3AC 


where  A C  *  C  -  C  ,  Ar  =*  r  -  r  and  the  signum  function  is  defined  as 
n  n  t  t  t 

n  n 

follows; 


1  Ac  >  0 

n 

sgn(ACn)  =0  ACn  =  0 


-1  AC  <  0 

n 


It  is  known  that  the  average  scattered  field  is  specular  [9,8]  .  This 

means  that  after  averaging  .t  (r)  over  the  stochastic  height  C  .  the  re- 

q 

suiting  dependence  of  the  averaged  current  on  r  is  of  the  form 

-►  *►  A  A 

exp(-  j  k .  *r  )  where  k.  =  kj  x  +  k.  y  is  the  transverse  part  of  the  inci- 
v  J  i  t'  i  ix  iy  v 

dent  wavenumber  vector.  In  general,  one  cannot  infer  this  particular  func¬ 
tional  dependence  before  averaging  because  C  is  a  stochastic  variable,  i.e. 

C(r  )  .  However,  in  the  case  of  (8)  where  C  and  C  ,n=],2,...,N  ,  are 
t  n 

taken  to  be  Independent  of  rf  (because  (8)  represents  a  discrete  approxima¬ 
tion  to  the  continuous  surface)  the  following  form  is  correct; 


J  (r)  *  j  (C)  exp  (-  j  k  •  rt> 

4  t 


(10) 


Substituting  (9)  and  (10)  in  (8)  and  performing  the  integration  over  r 
yields  the  following  result; 
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N 

j  (C)exp  (-jkit*rt)  =  Jq*  exp(-jki  •rt~jk1  Q  -  j  (^)  sgn  (Ar^ ) 

t  2  n=l 


*  J  exp{"  i  Vko  "  i 2  !A^J  "  j^'^L  +  Mkx-kit)-^t  Jdkx 

(11) 

i  A  A  a  *. 

where  J  =  2H  q  •  [  z  *  h)  ,  S  (• )  is  the  Fourier  transform  of  the  ntj*  support, 

qo  o  1  n 

and  r.-  is  the  vector  distance  from  the  origin  to  the  centroid  of  the  nc— 
tf‘o 

flat  area  projected  on  the  z = 0  plane.  The  curvatures  and  slopes  of  the  continu¬ 
ous  surface  are  taken  to  be  so  small  that  each  of  the  discrete  flat  areas 

A  ,n*l,2, •• *,N,  have  linear  dimensions  which  are  large  in  terms  of  a  wave¬ 

length.  Thus,  the  Fourier  transforms  of  the  support  functions  (S^)  will  be 
very  peaked  about  kA =■  k^  and  they  are  approximated  by  delta  functions.  Accom- 
plishing  the  kx-integration  and  dropping  the  common  exp  (-  j  k^*r^)  dependence 
yields 

N 

J(C)  ~  Jq*  exp(-jkizO  -  ^  j(Cn)sgn(A^n)  exp  (-j  kiz  |ACn|)  (12) 

n=l 


which  is  the  equation  that  must  be  solved  to  find  the  current. 

Before  continuing  on  with  the  solution  of  (12),  it  is  wise  to  review  the 

essential  rationale  which  permitted  the  simplification  of  (6)  to  (12).  The 

first  step  entails  the  observation  that  the  current  at  the  point  (x,y,0  may 

be  effectively  determined  by  its  interaction  with  the  current  at  an  unspecified 

number  of  other  points  (x  ,y  ,£  ),  n=l,2,**-,N,  on  the  surface.  In  view  of  the 

n  n  n 

very  slowly  undulating  nature  of  the  assumed  surface,  these  N  interacting  points 

are  approximated  by  large  flat  areas  (see  Figure  1).  The  elevations  of  these 

areas,  namely  £  and  5n,n-l,2 , • • • ,N,  are  random  but  independent  of  the  transverse 

spatial  coordinates  r£  and  r(  .  As  a  result  of  this  and  the  specular  nature 

n 

of  the  average  scattered  field,  the  dependence  of  the  current  upon  the  transverse 
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coordinates  r  may  be  inferred  directly  as  in  (10).  Finally,  the  height  depen¬ 
dent  part  of  the  current  is  found  to  be  a  solution  of  (12).  In  essence,  (12) 
relates  the  current  on  the  area  at  an  elevation  £,  to  the  currents  on  the  areas 
at  elevations  C  ,n=l,2, • • • ,N. 

As  a  first  step  toward  solving  (12)  it  should  be  noted  that  if  there  are 
no  additional  areas  of  interaction,  i.e.  N  =  0,  or  if  tne  other  areas  are  at  the 
same  elevation  as  the  area  at  C,  then 

j(0  =  Jqo  exp(- j  kiz  O  (13) 

which  is  the  correct  result  for  a  single  randomly  elevated  plane.  Of  particu¬ 
lar  significance  here  is  the  fact  that  (13)  is  independent  of  exactly  how  the 
limit  of  a  single  plane  is  achieved.  That  is,  the  proper  limit  is  obtained 
either  by  allowing  no  additional  areas  (N  = 0)  or  by  forcing  all  of  the  inter¬ 
acting  areas  to  coalesce  into  one  plane  (c  ■+■  ^,n=l,2 ,  •  •  •  ,N)  .  Equation  (13)  is 
also  recognized  to  be  the  physical  optics  or  single  scatter  approximation  for 
the  height  dependent  part  of  the  current. 

The  most  straightforward  approach  to  solving  (12)  is  to  start  with  one 
additional  area  (N  =  l)  and  then  build  on  this  result.  For  N  =  l,  (12)  becomes 

j  (O  =  exp  (-j  kiz  C)  -  j  C^)  sgn  (C  -  ^  exp  (-j  kiz  !?  -  ^  ! )  (14) 

which  relates  the  current  on  the  area  at  z  = £  to  the  current  induced  by  the 
incident  field  and  the  ci rrent  on  the  area  at  z  =  •  Rewriting  (14)  for  Liu- 

two  regions  C,  -  >  0  and  t,  -  <  0  yields 

j(?)  “  Jqg  exp(-j  k±z  C)  -  jC^)  exp  (-j  kiz  [t  -  ;11)  .',-^>0  (14a) 

and 

j(0  “  Jq^  exp(-j  kiz  c)  +  j(51)  exp  (j  kiz  -r^l)  r. -.'^<0  (14b) 

By  inspection  of  .ose  two  equations,  the  following  solution  is  obtained; 
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(15) 


j(0  -  Jq*  exp(- jklz  OlH-ACj) 

where  A£^  =  C  -  and  U(*)  is  a  unit  step  function  defined  as  follows 

0  A£  >  0 

UC-Aq)  =  1  (16) 

i  acl  <  0 

(the  fact  that  U(- Ac^)  =  1  for  Ac^ = 0  follows  from  the  appearance  of  the 
signum  function  in  (14)  and  its  definition  at  AC^=0).  In  verifying  that 
(15)  does  indeed  satisfy  (14),  one  must  be  careful  to  reverse  the  roles  of  C 
and  when  substituting  for  j(C^)  ,  i.e. 

J(Sl>  *  Jqo  e*P(-  jki2VU  (A^) 

Having  found  the  solution  of  (12)  for  N=1  ,  it  is  a  straightforward  task 
to  verify  that  adding  more  interacting  areas  simply  results  in  multiplying  (15) 
by  more  unit  step  functions.  That  is,  the  solution  of  (12)  is  as  follows; 

j(C)  -  J  q  exp(-jkizC)  TT  u(Cn  -  O  (17) 

n=*l 

This  is  an  interesting  and  not  altogether  unexpected  result.  According  to 
(17)  only  the  lowest  elevated  area  will  support  a  nonzero  current.  For 
example,  if  the  ^-boundary  is  the  lowest  then  A?n  >  0  for  n=l,2,...,N  and 
j(C)  is  just  the  current  due  to  the  incident  field.  Conversely,  if  is 

the  lowest  then  j(C^)  will  equal  the  current  due  to  the  incident  field  and 
all  of  the  other  boundaries  will  support  zero  current.  If  there  is  no  current 
on  a  boundary  then  the  boundary  produces  no  scattered  field  and  so,  for  all 
intents  and  purposes,  the  boundary  does  not  exist.  Thus,  (17)  predicts  that 
for  any  realization  of  the  random  boundaries,  only  one  (the  lowest)  deter¬ 
mines  the  scattering  from  the  surface.  This  result  is  entirely  reasonable 
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because  in  discretizing  the  original  integral  equation,  i.e.  (6),  Lhe  problem 
was  recast  into  N  +  l  interacting  perfectly-conducting  large  areas,  each  having 


all  the  characteristics  of  a  conducting  half-space  boundary.  Having  accom¬ 
plished  the  mathematical  solution  correctly  it  is  found  that,  for  all  intents 
and  purposes,  only  one  half-space  boundary  exists  for  any  realization  of  the 
randomly  elevated  surface.  Thus,  the  only  consequence  of  (17)  which  does  not 
have  an  immediate  physical  explanation  is  the  fact  that  it  is  always  the 
lowest  boundary  which  gives  rise  to  the  scattered  field.  Unfortunately,  no 
physical  explanation  for  this  particular  result  has  been  found. 

The  Average  Scattered  Field 

To  the  degree  of  approximation  provided  by  the  Fraunhoffer  diffraction 

•+ 

integral,  the  far-zone  scattered  field  in  the  direction  kg  and  at  a  distance 
R  from  the  mean  surface  is  given  by 

JC?t  )exp(jks  -Tt  +  jkszCm)d^t  i  (18) 
m  t  m  "  m 


where  k  =k  +  k  z  *  k  k  ,  n  is  the  impedence  of  free-space  and  g  is  given 

3  St  SZ  O  S  O 

by  (3).  The  sum  is  over  all  M+l  flat  areas  comprising  the  total  discretized 

surface  and  (£0>A0)  represents  the  flat  area  having  a  height  C  and  area  A. 

2 

Consistent  with  the  previous  stipulation  that  A  »  A  ,  (18)  reduces  to  the 

TO  O 

following 


A  A  A  A 


E  -  j  2trk  n  H  g(R)  [k  xk  x(z*h)]6(k  -k.  )  6  (k  -k.  ) 

s  J  o  o  o  °  s  s  sx  IX  sy  ly 


M 


j(CTn)expfj(ksz-kiz)CmJ 


(19) 


If  £  ■  0,m=0,l,2,  ••• ,  this  expression  for  E  is  too  large  by  the  factor  M  +  l; 
in  s 

this  error  results  from  treating  each  Am  as  essentially  infinite  and  it  can  be 
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rectified  by  dividing  the  rhs  of  (19)  by  M  +  l.  Taking  the  average  of  the  cor¬ 
rected  version  of  (19)  and  noting  that  the  average  of  j (C  )exp(ik  C  )  is 

id  sz  m 

independent  of  m  yields 


<E>*j2iTk  n  K  g(R)  [k  *  k  *(zxh)]6(k  -k.  )  6  (k  -k .  ) 
s  J  ooo°  s  s  sx  ix  sy  iy 


•<J(0  exp[j(ksz-kiz)  C  ]> 


(20) 


where 


exp(-jklzC)  N  =  0 

J(0  - 

N 

exp  (-jk  C)  IT  U(5  -O  N  >  1 

iz  n*i  n  — 


(21) 


The  factor  multiplying  the  <•>  term  in  (20)  is  identical  to  the  average  scat¬ 
tered  field,  Eg°,  for  the  case  of  a  plane  located  at  z  =  0  in  the  far-field 

approximation.  Thus,  <E  >  has  the  same  direction  as  E  and  it  differs  in  com- 

s  s 

plex  amplitude  by  the  term  <•>  in  (20).  Since  the  scattered  field  is  specular, 

k  ■  ±  k ,  and  for  the  field  in  the  upper  half-space  (free-space)  k  *  -  k  ; 

SZ  1  z  sz  xz 

thus,  (20)  becomes 

<V7V  =  <0  «>  e*P  l  -  j  2kiz  (22) 


where  k,  =  -k  cos  0  . 
iz  o 


For  N»0,  (22)  yields  the  following  result; 


(23) 


where  p^  is  the  characteristic  function  of  C  .  For  N  >  1  ,  there  results 
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<E  /E°> 
s  s 


00 

r 

_  00  OO  to  — \ 

r  r  r 

1  exp (j  2|klg|  O  < 

J 

I  ...  Pte1,r2,...,ru,r)dcldr2...drN> 

-OO 

1  J  J  J  \ 

-c  c  c 

dr 


(24) 


where  P(^i»^2’ '  *  *  t*le  ^0^nt  probability  density  function  of  the 

indicated  random  variables.  This  result  for  the  relative  strength  of  the  aver¬ 
age  field  involves  two  parameters  -  the  degree  of  correlation  between  the 
non-overlapping  areas  and  the  number  of  areas  which  are  important  to  the  scat¬ 
tering.  According  to  the  similarity  obtained  in  [8}  and  reviewed  in  the  first 
section  of  this  paper,  total  correlation  represents  the  gently  undulating  sur¬ 
face  (£  >>  A  )  while  complete  decorrelation  now  represents  the  uniformly  rough 
surface  (£«Aq)  limit.  When  the  areas  are  all  correlated,  it  can  be  shown 
using  conditional  densities  that  (24)  reduces  to  (23)  independent  of  N.  This 
is  as  it  should  be  because  for  £  »  Aq  (24)  represents  the  case  of  essentially 
a  randomly  elevated  plane.  When  the  non-overlapping  areas  are  uncorrelated, 

(24)  is  very  highly  dependent  upon  N  as  should  be  expected  since  this  repre¬ 
sents  the  case  of  very  strong  mutual  interaction  on  the  surface.  This,  of 
course,  is  simply  a  reiteration  of  the  fact  that  as  the  flat  areas  become 
decorrelated  the  mutual  interaction  effects  on  the  surface  increase.  The  re¬ 
mainder  of  this  paper  will  be  devoted  to  understanding  and  quantifying  the 
importance  of  N  for  the  case  of  uncorrelat*.  1  surfaces  (£<<Aq). 

For  an  arbitrary  joint  density  function  representing  uncorrelated  heights, 
it  is  difficult  to  simplify  (24).  What  is  happening  here  is  that  the  mathe¬ 
matics  of  probability  theory  are  beginning  to  cloud  the  physics  of  the  process. 
To  overcome  this  limitation,  consider  the  case  where  decorrelation  implies 
statistical  independence.  For  a  statistically  independent  process,  the  joint 
density  factors  into  a  product  of  the  marginal  densities,  i.e. 

p(c1,c2,...,cn.o  -  p(;1)p(c2).  ..pUN)pU) 
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The  Cn- integration  in  (24)  may  be  rewritten  as  follows; 

C 

1  -  F(C) 


r  ? 

J  '  Vdt»  ■ 1  -  /  <>  V5, 


where  F(£)  is  the  common  distribution  function  for  C,C, ,C  ,  ;  that 

1  2  N 

is,  all  the  random  elevations  are  identically  distributed  because  of  the 
statistical  homogeneity  stipulation.  Thus,  when  decorrelation  implies  sta¬ 
tistical  independence,  (24)  simplifies  to  the  following  form; 


<E  /E°> 
s  s 


exp(j  2|klz|  C)P(OU-F(?)]MdC 


(25) 


Since  F(-°°)  *  0  and  F(+°°)  «  1  ,  the  effect  of  the  N  interacting  areas  on 
the  surface  is  to  skew  the  integrand  toward  the  negative  range  of  £  in  (25). 
Unfortunately,  it  is  difficult  to  proceed  beyond  this  general  result  without 
a  specific  form  for  p(Q  and  F(£)  for  two  reasons.  First,  it  is  diffi¬ 
cult  if  not  impossible  to  find  a  meaningful  form  for  the  integration  in  (25)  for 
arbitrary  p(0  and  F(£).  Second,  and  more  importantly,  is  the  simple  fact  that 
there  is  no  formula  for  determining  the  effective  number  of  interacting  areas  N; 
it  must  be  determined  by  computing  (25)  as  a  function  of  N  and  rationalizing 
the  results  with  the  physics  of  the  rough  surface  scattering  process. 


Numerical  Results  and  Interpretation 

For  Gaussian  surfaces  (the  only  class  of  surface  statistics  for  which  the 
similitude  has  been  proven  [8]), the  marginal  density  is 

p(0  »  - 1 - exp(-  £2/2<£2>) 

V2tt<C2> 

and 
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(26) 


1  -  F(0  -  \  erfc^//2<C2>  ) 

2 

where  <£  >  is  the  mean  square  surface  height  and  erfc(*)  is  the  complemen¬ 
tary  error  function.  Substituting  these  forms  in  (25)  and  introducing  the 
normalized  variable  n  ■  c/*4<£;2>  yields; 


<E  /E°>  = 
s  s 


1 

/a 


CO 

I 


exP(j  2/F  E  n  -  n2)  [erfc(n)/2]N  dn 


(27) 


where  I  *  jk^  \Y<t?>  -  k0*^C2>  cos  0  is  the  Rayleigh  parameter.  For  arbi¬ 
trary  E  ,  this  integration  is  straightforward  but  must  be  accomplished 
numerically.  For  E*0  ,  the  integral  can  be  shown  to  yield  1/(N+1); 
consequently,  as  N  increases  (27)  monotonically  decreases.  Figure  2 
illustrates  the  behavior  of  the  magnitude  of  the  average  scattered  ield  as  a 
function  of  E  for  Na 0,1,2,  and  3-  Of  particular  note  here  is  the  fact 
that  the  curves  begin  to  flatten  out  as  both  E  and  N  increase.  This  figure 
also  shows  a  very  interesting  result;  for  N-0,  1,  and  2  there  is  a  specific 
range  of  E  values  for  which  each  curve  provides  the  maximum  average  scattered 
field  compared  to  all  the  other  curves  (for  N  =  0  the  range  is  0  £E  ~  1  ).  Fig¬ 
ure  3  shows  a  plot  of  N  values  as  a  function  of  the  corresponding  E  range 
over  which  the  average  scattered  field  is  maximized.  That  is,  for  each  value 
of  N  and  E  corresponding  to  a  point  in  Figure  3,  equation  (27)  is  maximized, 
i.e.  any  other  N  will  provide  a  smaller  value  for  (27). 

In  order  to  appreciate  the  implication  of  this  result  it  is  necessary  to 
recall  the  meaning  of  N.  As  simulated  in  this  analysis,  N  represents  the  effec¬ 
tive  number  of  uncorrelated  areas  on  the  surface  which  interact  with  the  refer¬ 
ence  point  on  the  surface.  Equation  (27)  describes  the  influence  of  these 


N-area  interactions  on  the  coherent  scattered  field  as  a  function  of  the 


Rayleigh  parameter  or  surface  roughness.  From  a  physical  point  of  view,  one 
should  expect  that  for  small  roughness  each  point  on  the  surface  scatters 
independently  of  all  the  other  points;  hence,  this  corresponds  to  the  N  =  0  or 
no  interaction.  As  the  surface  becomes  rougher,  it  is  certainly  reasonable 
to  expect  that  the  region  of  interaction  on  the  surface  should  increase;  hence, 
N  should  also  increase.  With  this  bit  of  simple  physics  in  mind,  it  is  now 
possible  to  interpret  the  results  in  Figures  2  and  3.  In  essence.  Figures  2 
and  3  show  the  range  of  Z  values  for  which  a  given  number  of  interacting  areas 
on  the  surface  are  dominant  or  are  most  important  to  the  scattering  process. 
Figure  3  shows  that  for  £~5  ,  N  varies  very  nearly  as  £  .  Hence,  the  number 
of  interacting  points  or  regions  on  the  surface  is  increasing  as  the  cube  of 
the  roughness.  To  the  author's  knowledge,  this  is  the  first  time  that  such 
a  direct  correspondence  has  been  demonstrated. 

It  should  be  noted  that  the  parameter  N  may  be  interpreted  [11]  somewhat 
differently  than  presented  above.  Although  this  secondary  interpretation 
does  not  follow  directly  from  the  analysis  presented  here,  it  would  certainly 
seem  to  apply  to  the  original  uniformly  rough  surface.  In  particular,  the 
parameter  N  may  be  thought  of  as  the  number  of  additional  reflections  that  an 
incident  ray  undergoes  due  to  multiple  scattering  on  the  surface.  From  this 
point  of  view,  each  value  of  N  is  clearly  only  dominant  over  a  limited  range 
of  roughness  [11].  Furthermore,  this  analogy  provides  an  alternate  reason¬ 
ing  for  moving  from  the  N=nQ  to  the  N*nQ+l  curve  in  Figure  2  as  the  Rayleigh 
parameter  increases.  That  is,  once  the  N*no  curve  drops  below  the  N  =  no  + 1 
curve  in  Figure  2,  this  signifies  the  need  to  include  more  multiple  scattering 
or  ray  reflections  on  the  surface.  This  interpretation  and  the  results  in 
Figure  3  imply  that  the  order  of  multiple  scattering  or  the  number  of  ray 
reflections  on  the  surface  increases  as  the  cube  of  Rayleigh  parameter. 

Clearly,  this  interpretation  is  synonymous  with  the  increasing  area  of 
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interaction  on  the  surface. 

Given  the  above  rationale  for  selecting  N  as  a  function  of  Z,  how  then 
does  one  construct  a  continuous  curve  of  the  average  scattered  field  as  a 
function  of  Z  ?  A  technique  that  has  been  found  to  be  relatively  fast  is  the 
following.  One  chooses  three  contiguous  values  of  N  and  computes  (27)  as  a 
function  of  Z  for  each  of  the  three  values.  One  then  finds  the  value  of  Z 
which  maximizes  (27)  for  the  middle  N  and  this  provides  one  point  on  the 
plot  of  average  scattered  field  strength  versus  Z  .  This  process  is  continued 
until  a  smooth  curve  can  be  drawn  through  the  points.  Because  N  varies  as 
the  third  power  of  I  ,  the  maxima  of  (27)  are  very  close  together  and  so  not 
many  points  are  required  before  the  trend  of  the  maxima  can  be  established. 

The  only  point  where  there  will  be  some  degree  of  interpolation  required  is 
near  the  transition  from  N-0  to  N»l. 

Figure  4  compares  (27),  as  computed  according  to  the  above  prescription, 
with  the  results  of  DeSanto  and  Shisha  [6].  (Both  results  clearly  show  a 
larger  average  scattered  power  than  is  predicted  by  physical  optics.)  The 
disagreement  is  due  most  probably  to  the  fact  the  integral  equation  for  the 
average  scattered  field  used  by  DeSanto  and  Shisha  is  only  the  first  term  in 
an  infinite  hierarchy  of  integral  equations.  The  solution  represented  by  (27) 
is  also  an  approximation  because  the  large  flat  areas  comprising  the  substitute 
surface  have  essentially  been  replaced  by  planes.  However,  this  simplifica¬ 
tion  does  not  significantly  alter  the  basic  processes  responsible  for  the 
effects  of  multiple  interaction  on  the  average  scattered  field  and  so  it  is 
reasonable  to  express  a  high  degree  of  confidence  in  these  results. 

Conclusions 

The  primary  purpose  of  this  paper  is  to  present  an  alternative  approach 
to  the  problem  of  coherent  scattering  from  a  uniformly  rough  surfaces,  i.e. 
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surfaces  for  which  the  correlation  length  is  significantly  smaller  than  the 
electromagnetic  wavelength.  The  motivation  for  this  effort  is  the  need  to 
obtain  a  degree  of  physical  insight  which  is  not  present  in  the  exact  singular 
integral  equation  describing  this  problem.  The  approach  is  based  upon  the 
remarkable  similarity  [8],  with  respect  to  the  intermediate  mathematical  re¬ 
sults,  between  the  stochastic  Fourier  transform  of  the  total  current  on  the 
gently  undulating  surface  and  the  part  of  the  stochastic  Fourier  transform  of 
the  current  on  the  uniformly  rough  surface  which  is  responsible  for  the  aver¬ 
age  scattered  field.  This  similarity  leads  to  the  replacement  of  the  uniformly 
rough  surface  by  a  discretized  version  of  the  gently  undulating  surface  com¬ 
prising  large,  flat,  non-overlapping,  un correlated  areas.  This  substitution 
is  justified  by  the  fact  that  the  stochastic  Fourier  transform  of  the  current 
on  the  discretized  surface  obeys  exactly  the  same  governing  relationships  as 
for  the  uniformly  rough  surface. 

A  second  major  step  in  the  analysis  entails  identifying  the  effective 
number  of  interacting  areas  on  the  discretized  surface  in  terms  of  orders  of 
multiple  scattering  on  the  surface  and  recognizing  the  necessity  to  monotoni- 
cally  increase  the  order  of  multiple  scattering  with  increasing  surface  rough¬ 
ness.  In  particular,  it  is  found  that  the  number  of  interacting  areas  on  the 
surface,  or  alternatively  the  order  of  multiple  scattering,  necessary  to  maxi¬ 
mize  the  average  scattered  field,  varies  as  the  cube  of  the  Rayleigh  rough¬ 
ness  parameter.  The  maximum  property  of  the  average  scattered  field  is  con¬ 
sistent  with  earlier  arguments  [8]  that  the  uniformly  rough  surface  produces 
a  maximum  coherent  field  for  a  fixed  Rayleigh  parameter. 

The  essence  of  the  approach  is  that  of  solving  a  simpler  problem  which 
can  be  shown  to  have  the  same  mathematical  properties  with  respect  to  the 
average  scattered  field  as  the  actual  surface.  Clearly,  it  is  desirable 
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t  compare  the  results  obtained  herein  with  a  solution  of  the  exact  Integral 
equation  for  the  average  scattered  field.  However,  this  must  await  further 
study  of  the  integral  equation  and,  in  particular,  the  analytic  properties 
of  the  average  scattered  field  for  "off-shell"  (k  ^-k.  )  conditions. 
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Figure  la.  Geometry  of  the  rough  surface  interface. 

Figure  lb.  Discrete  approximation  to  the  true  rough  surface  in  Figure  la. 

The  dashed  continuous  line  is  the  original  surface  while  the  non¬ 
overlapping  horizontal  line  segmets  form  the  discrete  approxima¬ 
tion  to  the  surface. 

Figure  2.  A  plot  of  the  coherent  scattered  power  as  a  function  of  surface 

roughness  for  N  =  0,l,2,&3.  Note  that  for  each  N,  there  is  a  range 
of  surface  roughness  for  which  the  coherent  power  is  larger  than 
for  all  other  values  of  N. 

Figure  3.  A  plot  of  the  N  value  which  maximizes  the  average  scattered  power 
over  the  indicated  range  of  surface  roughness.  The  figure  illus¬ 
trates  the  increasing  region  of  multiple  interaction  on  the  surface 

as  the  surface  roughness  increases.  The  solid  line  represents  the 
3 

relationship  N  *  Z  . 

Figure  4.  A  comparison  of  the  coherent  scattered  power  as  a  function  of  sur¬ 
face  roughness  as  predicted  by  single  scattering  (N  =  0),  DeSanto's 
&  Shisha’s  computations,  and  equation  (27).  N,  or  the  number  of 
interacting  points  or  regions  on  the  surface,  has  been  replaced  by 
a  continuous  function  of  £  which  fits  the  results  in  Figure  3. 
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Abstract 

This  paper  develops  new  results  pertaining  to  electromagnetic  scatter¬ 
ing  from  perfectly  conducting  random  surfaces  for  which  decorrelation  Joes 
not  imply  statistical  independence.  Using  an  exact  theory  for  Lhe  current 
induced  on  the  surface  and  the  far-field  approximation  for  the  scattered 
field,  it  is  shown  that  the  incoherent  scattered  power  consists  of  two  parts. 
The  first  part  corresponds  to  the  so-called  diffuse  scattered  power.  The 
second  part  is  specular  in  its  angular  dependence  and  is  a  direct  consequence 
of  the  fact  that  the  two  point  joint  density  for  the  surface  height,  slopes, 
etc.  does  not  reduce  to  the  product  of  the  single  point  joint  densities  for 
infinite  separation  distances.  Computations  for  a  gently  undulating  jointly 
exponentially  distributed  surface  show  that  the  incoherent  specular  power  is 
equal  to  coherent  scattered  power  for  the  Rayleigh  parameter  near  one.  When 
the  Rayleigh  parameter  is  large,  the  incoherent  specular  power  is  signifi¬ 
cantly  larger  than  the  coherent  power.  The  analysis  further  indicates  that 
scattering  measurements  provide  an  ideal  way  for  identifying  this  class  of 
surfaces. 

Introduction  and  Summary 

Frequently,  when  comparing  measurements  of  scattering  from  randomly 
rough  surfaces  with  theoretical  models  there  are  little  or  no  data  on  the 
statistical  properties  of  the  surfaces.  This  limitation  means  that  the 


surface  statistics  must  be  assumed.  Unfortunately,  t he  choice  ot  the  .surface 
statistics  is  usually  based  upon  analytical  convenience  rather  than  reality. 

If  the  assumed  statistics  are  not  a  good  approximat  ion  to  the  actual  siiil.ue, 
this  can  lead  to  erroneous  rejection  of  the  theory  or  expensive  additional 
experimentation. 

Recently,  Lennon  and  Papa  [1980,  1981]  employed  digitized  terrain  maps 
to  obtain  estimates  of  very  large  scale  surface  height  sLat isL ics  which  could 
then  be  used  to  predict  scattering  from  the  region  in  question  [Papa,  e_t  al. , 
1980].  Because  of  the  complexity  of  the  theoretical  scattering  model,  it  is 
not  practical  to  use  measured  statistics  such  as  joint  height  histograms 
directly  in  the  scattering  computations.  A  reasonable  alternative  is  to  "l>cst- 
fit"  these  statistics  to  a  functional  form  for  the  joint  height  probability 
density  function  which,  in  turn,  is  amenable  to  use  in  the  scattering  theory. 
Since  the  functional  forms  for  the  joint  height  statistics  may  be  quite  gener¬ 
al,  scattering  theories  based  upon  Gaussian  surface  statistics  are  no  longer 
adequate.  Beckmann  [1973]  has  considered  the  problem  of  scattering  from  non- 
Gauasian  surfaces  when  the  only  height  data  available  are  the  marginal  den¬ 
sity  and  the  correlation  function.  However,  his  analyses  and  results  are 
restricted  to  the  case  where  decorrelation  of  the  surface  height  implies  sta¬ 
tistical  independence. 

The  purpose  of  this  paper  is  to  consider  scattering  from  surfaces  lor 
which  decorrelation  does  not  imply  statistical  independence.  The  motivation 
for  this  study  is  threefold.  First,  there  is  no  obvious  reason  why  a  real 
surface  should  be  statistically  Independent  when  it  is  decorrelaied.  Also, 
from  a  mathematical  point  of  view  statistical  independence  is  not,  in  general, 
implied  by  a  lack  of  correlation  [Papoulis,  1965].  Finally,  there  are  a 
relatively  large  number  of  joint  probability  functional  forms  for  which 
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decorrelation  and  statistical  independence  are  not  synonymous  [Lennon  and 
Papa,  1981]  and  there  is  no  apparent  reason  why  these  should  be  rejected 
a  priori  when  attempting  to  determine  the  "best  fit"  functional  form. 

The  analysis  starts  with  describing  the  far-zonc  scattered  field  by  tiie 
standard  far-fleld  diffraction  integral  Involving  the  current  Induced  on  the 
surface  by  the  incident  plane  wave  and  the  Fourier  phase  kernel.  Averages 
necessary  to  compute  the  mean  scattered  field  and  its  variance  are  converted 
to  convolutions  in  a  non-stochastic  Fourier  transform  space  [Brown,  1982]. 

The  advantage  of  this  approach  is  that  the  spatial  dependence  of  the  trans¬ 
formed  current  is  known  and  this  in  turn,  implies  that  the  average  scattered 
field  is  specular,  e.g.  it  is  nonzero  for  one  and  only  one  scattering  angle. 
The  second  moment  of  the  scattered  field  contains  a  specular  component  due  to 
the  nonconvergent  nature  of  the  two  point  joint  density  for  the  surface 
height,  slopes,  curvatures,  etc.  of  the  surface  as  the  distance  between  the 
two  points  on  the  surface  approaches  infinity.  When  decorrelation  implies 
statistical  independence,  this  specular  component  is  exactly  canceled  by  the 
subtraction  of  the  square  of  the  mean  scattered  field;  Lhus,  the  variance  of 
the  scattered  field  contains  no  specular  components.  However,  if  the  surface 
is  not  statistically  independent  when  it  is  uncorrelated,  the  variance  of  the 
scattered  field  is  shown  to  have  a  nonzero  specular  component  which  is  depen¬ 
dent  upon  the  difference  between  the  decorrelated  joint  density  and  the 
square  of  the  marginal  or  single  point  density  for  the  surface  height,  slopes, 
curvatures,  etc.  This  result  is  a  simple  consequence  of  the  fact  that  no  two 
points  on  the  surface  scatter  statistically  independent  of  each  oilier  regard¬ 
less  of  the  distance  between  them. 

In  order  to  illustrate  this  effect,  a  very  gently  undutaiiug  surface 
which  is  amenable  to  the  physical  optics  approximation  for  the  Induced  surface 
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current  is  assumed.  The  two  polut  Joint  height  probability  density  i unction 


is  taken  to  be  exponential  and  the  specular  part  of  the  scattered  field  v.n  i- 
ance  is  compared  to  the  square  of  the  mean  scattered  field  as  a  function  of 
the  Rayleigh  roughness  parameter  (I).  The  specular  part  of  the  scattered 
field  variance  goes  to  zero  as  the  Rayleigh  approaches  zero;  however,  as  the 
Rayleigh  parameter  increases,  it  also  increases  to  a  peak  value  in  tin.'  neigh¬ 
borhood  of  E«1  where  it  is  equal  to  square  of  the  mean  scattered  field.  for 
larger  values  of  E,  the  specular  part  of  the  variance  greatly  exceeds  the 
square  of  the  mean  scattered  field;  the  former  behaving  asymptotically  as 

_  3 

E  while  the  latter  goes  as  E 

The  importance  of  this  analysis  rests  primarily  in  its  implications.  If 
the  Rayleigh  parameter  is  very  email  (E«l)  or  very  largo  (E  >>  1 ) ,  it  is 
doubtful  that  the  specular  component  of  the  variance  of  the  scattered  field 
will  be  measurable.  However,  if  the  geometric  and  electromagnetic  parameters 
of  the  experimental  scattering  system  can  be  adjusted  so  that  L  ~  1  then  this 
specular  part  of  scattered  field  variance  should  be  detectable.  If  no  specu¬ 
lar  part  of  the  field  variance  is  measured  then  one  may  safely  conclude  that 
the  surface  is  statistically  Independent  when  it  is  decorrelated.  A  final 
implication  of  this  analysis  is  that  one  should  exercise  great  care  in  tit- 
ting  arbitrary  functional  forms  for  joint  probability  densities  to  histogram 
data  for  purposes  of  predicting  theoretical  scattering  l rum  the  surface. 

Scattering  Analysis 

The  stochastic  surface  z  ■  £(x,y)  is  assumed  to  be  a  zero  mean  statis¬ 
tically  homogeneous  process.  Beneath  the  surface,  the  medium  is  taken  to  he 
perfectly  conducting  while  free-space  comprises  the  medium  above.  A  magnetic 

•+  a  4  a 

field  of  the  form  H  (r)  ■  H  h  exp(-j  k  *r)  is  incident  upon  the  surface;  h  is 

o  l 

A  -► 

the  polarization  of  the  magnetic  field,  k^  -  kj/k0  is  its  direction  of 


propagation,  and  kQ  *  2ir/AQ  is  the  wavenumber.  Although  the  far-field  scattering 
approximation  leads  to  difficulties  when  dealing  with  unbounded  illumination 
[Saucer,  1969],  these  problems  can  be  circumvented  by  relating  the  results  to 
flat  surface  "scattering"  in  the  same  approximation.  WiLh  this  caveat  in 
mind,  the  far-zone  scattered  electric  field  is  given  by  the  following  form; 
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where  nQ  ia  the  impedance  of  free-space,  g(RQ)  ■  exp(-  j kQR^)/Ru ,  Rq  is  the 
distance  from  the  origin  of  the  surface  based  coordinated  system  to  the  far- 
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field  point,  k  /k  *  k  is  the  direction  of  the  far-field  point  of  observa- 
r  s  o  s 

-f  -> 

tion,  and  r£  «  xx  +  yy  is  a  point  on  the  z*0  plane.  The  quantity  J(r)  is 


defined  as  follows; 


J(r)  -  J  (r)[l>Oc/dx)2  +  Oc/3y)2]' 


where  Ja<r)  is  the  current  density  on  the  surface  and  8t;/3x  and  9£/3y  are  the 
x  and  y  slopes  of  the  surface  at  the  point  (£,rt).  Redefinition  of  the  cur¬ 
rent  as  in  (2)  permits  an  integration  over  the  z»0  plane  in  (1)  rather  than 
the  actual  surface.  With  T  »  j  kQno  g(R0)/4n  ,  the  average  scattered  field 


is  given  by 


-  r  I  I" k  x < J(r)exp (j  k  ?,)>"]  exp(jk  *r  )dr 
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where  £  +  k  z  »  k  and  the  <•>  brackets  denote  an  ensemble  average  over 

8  SZ  S 

all  stochastic  quantities  upon  which  J(r)  depends.  With  p^Ct.Z)  equal  to 
the  single  point  joint  probability  density  function  for  the  stochastic  vari¬ 
ables  £  and  Z  ,  the  average  in  (3)  is  as  follows; 


<J (r)exp (jk  £) 
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J(r)exp(jkazOp1  (c ,Z)JcdZ 


(4) 


where  Z  is  symbolic  of  all  Lhe  surface  height  derivatives  wlii.lt  dote i mine 
J(r).  Since  this  number  is  in  general  infinite,  Z  stands  for  V 
[Brown,  1982]  and  the  dependence  of  these  higher  order  surface  derivatives 
upon  r  is  implied.  With 

■+  A  /\  ■+  'p  a  ■+  a  -► 

J (r)  -  k  *k  *J(r)  *  (k  *J)k  -J  (5) 
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and  using  (4),  the  squared  magnitude  of  the  average  scattered  field  is  easily 
shown  to  reduce  to  the  following  form; 
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•  expj^jk^  *(rt-rt' jjdCd^'dZdZ'd^dr^  (6) 

where  the  synbol  *  denotes  the  complex  conjugate  operation  and  p^(C',Z')  is 
not  conjugated  because  it  is  real.  Following  a  similar  development,  the 
mean  squared  scattered  field  may  be  expressed  as  follows; 


<  E. 


|r| 
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(rt-rt'  )1  d£dS'  dZdZ'dr^dr^ 


(7) 


where  Art  “  rt“rt'  •  The  function  P2  in  (7)  is  the  two  point  joint  probabil¬ 
ity  density  function  which  depends  upon  two  sets  of  stochastic  variables  and 
the  horizontal  vector  distance  between  them. 

If  the  stochastic  Fourier  transforms  are  defined  as  follows; 
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f(rt,k,K)  ~  J J 
(k,K)  -  JJ Pi<C, 


J(r)exp[jU  +  jKZjd^dZ 
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+  KZ  +  K'Z'JJdCdC'dZdZ' 


and  it  is  noted  that 


F[j*(rt,C,Z)]  -  ^J(tt,-k,-K)|*  (8d) 

where  F  symbolizes  the  stochastic  Fourier  transform  operation,  the  integra¬ 
tions  over  the  stochastic  variables  in  (6)  and  (7)  may  be  alternatively  ex¬ 
pressed  as  convolutions  in  transform  space  [Brown,  1982].  The  advantage  of 
this  approach  is  that  the  new  variables  of  integration,  i.e.  the  transform 
coordinates,  do  not  depend  upon  rt>  r  ,  or  Arfc  .  Rewriting  the  integrals 
In  (6)  and  (7)  as  noted  above  yields 


i<Vi2  - 


(2tt)2+2S“ 


Hi f t*k’K),[^t*“k,*‘K,)]  r,i(ksz"k*‘ 


<|eJ2>  - 


Pl(“kSZ~k.~K)}eXP[Jks  •(V^)]dkdk,dKllK'dV^ 


X2+2S„ 
(2ir)  00 


!>„(k  -k,-k  -k'.-K.-K'  ;Ar  ) 

2  SZ  sz  IJ 


exp^jk^  •  (r j.-r^jdkdk' dKdK' dr^dr^ 


4-7 


The  quantity  S  is  equal  to  £  i  and  it  results  from  the  dependence  of 

i-2 

3  (r)  on  all  orders  of  surface  height  derivatives  [Brown,  1982).  The  euher- 

0 

i  i2 

ent  scattered  power  is  proportional  to  |<Eg>|  while  ihc  incoherent 

.  ->  .  2  .  -*■  .  2 

scattered  power  is  proportional  to  Var(E  )=<[E|>-[^E>[  or  the 

s  s  s 

variance  of  the  scattered  field.  Using  (9)  and  (10),  the  variance  of  the 
scattered  field  may  be  written  as  follows; 
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In  view  of  (5)  and  (8a),  the  stochastic  Fourier  transform  of  j(r)  is 
proportional  to  the  stochastic  Fourier  transform  of  J(r).  However,  it  has 
been  shown  [Brown,  1982]  that,  in  order  to  satisfy  the  stochastic  Fourier 
transformed  magnetic  field  Integral  equation  for  the  current,  J(rt>k,K)  must 
have  the  following  dependence  upon  rt; 


3 (r ^ , k , K)  -  j(k,K)exp(-jki^*rt) 


(12) 


where  k,  is  defined  by  k,  »  k.  +  k.  z  .  That  is,  the  stochastic  Fourier 
it  i  it  iz 

transform  of  the  slope  normalized  current  exhibits  the  same  dependence  upon 
rfc  as  the  Incident  magnetic  field.  Using  (5),  (8a),  and  (12)  in  (11)  and 
converting  to  a  Arfc  coordinate  yields 


Var(Es)  - 


“Sr  Iff  If 


..  .  2+2S 
(2ir)  * 


*{p2<k8z"k*_ksz"k,»"K’~K,;A^t)  'Pl(ksz"k'“K)  ~VL  (  ks/.~k*'K)3 
•  exp£jArt*(ks  dkdk'dKdK'dArt 


4-8 


where  A  is  the  illuminated  area  and  the  summation  is  over  the  values  q-x,y, 
and  z. 

It  should  be  noted  that  (13)  is  an  exact  result  Involving  no  approxima¬ 
tions  other  than  the  far  field  integral  description  of  the  scattering  pro¬ 
cess.  The  function  is  also  called  the  single  point  joint  character- 

istic  function  of  the  stochastic  variables  ?  and  Z  while  p-,  (  4 ,  •  , 4  , 4  ;  Ar  )  is 
the  two  point  joint  characteristic  function  of  C,£',Z,  and  Z' .  The  important 
point  to  note  about  (13)  is  that  apart  from  the  phase  factor  only  p7  depends 
upon  Ar^ .  The  explicit  dependence  of  p^  on  Ar^  comes  abouC  through  the  sur- 
face  height  correlation  fvmction  R(Ar^)  and  its  higher  order  derivatives  where 

R(A^t)  -  <C(rt' +A^t);(^t')> 

When  |Art|-*-°°  ,  the  correlation  function  and  its  higher  order  derivatives  ap¬ 
proach  zero  and  this  corresponds  to  the  condition  known  as  decorrelation.  If 
decorrelation  implies  statistical  Independence  then 

lim  P2  "  Pi^-k.-K)?  (-k  -k',-K')  (14) 

|Art|-K» 


or  the  two  point  characteristic  function  factors  into  a  product  of  single 

point  characteristic  functions.  When  this  occurs,  the  Art-iutegrand  in  (13) 

vanishes  as  | Arfc |  and  so  it  has  a  finite  support.  This  implies  that  the 

Var(E  )  will  contain  no  delta  function  dependency  upon  k  -k  .  Thus,  when 
8  st  H 

decorrelation  implies  statistical  independence,  Var(E  )  will  contain  no  spec- 

b 

ular  terms  because  specularity  in  the  far-field  approximation  implies  a  terra 
.  +  + 

which  varies  as  o(k  -k  )  . 

at  it 

If  decorrelation  does  not  imply  statistical  independence  then  (14)  is 


Invalid  and  the  term 


will  be  nonzero.  The  consequences  of  this  result  may  !»«.•  i-.isily  sron  l>y  ivgu- 
larlzing  the  Art~integrand  in  (13)  by  subtracting  and  then  adding  terms  con¬ 
taining  the  above  limiting  form.  This  operation  leads  to  the  following 
result; 
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(2ir) 
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[1-k  -q]1  (k,K)j*(-k\-K’) 
o  q  q 


*  ^P2  (k8r“k* -k3z-k'.-K,  -K';  I  Art  |  «»)-p1  (kgz-k  ,-K)^  (-ksz-k ' ,  -K ’  )^  dkdk  ’ dKdK  ’ 

(15) 

The  second  term  in  (15)  is  clearly  specular  in  its  angular  dependence  and  it 
is  entirely  a  consequence  of  the  difference  between  decorrelation  and  statis¬ 
tical  independence. 

It  might  appear  that  this  result  is  a  freak  consequence  of  the  unbounded 
nature  of  the  incident  field;  however,  such  is  not  the  case.  The  analysis  may 
be  carried  through  using  a  plane  wave  having  a  finite  support  (a  beam  plane 
wave)  provided  that  the  linear  dimensions  of  the  illuminated  area  are  large 
relative  to  the  decorrelation  length  of  the  Burface.  The  only  essential  dif¬ 
ference  between  the  limited  area  result  and  (15)  will  be  the  appearance  of 
the  diffraction  pattern  of  the  finite  illuminated  area.  If  the  linear 
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dimensions  of  the  illuminated  area  are  not  large  compared  to  the  decorrelation 


length,  the  above  results  do  not  apply.  However,  this  is  an  entirely  dif¬ 
ferent  problem  which  will  not  be  considered  here. 

Discussion 

The  appearance  of  a  specular  term  in  the  expression  for  the  incoherent 
power  is  certainly  a  bit  surprising.  Consequently,  it  is  desirable,  if  not 
mandatory,  to  seek  some  physical  explanation  for  this  result.  Unfortunately, 
no  such  explanation  has  been  found  primarily  because  the  statistical  moments 
of  the  scattered  field  are  obtained  by  mathematical  operations  governed  by 
the  laws  of  probability  rather  than  physics.  From  a  very  fundamental  point 
of  view,  this  problem  is  closely  akin  to  inquiring  into  the  physical  signifi¬ 
cance  of  statistical  dependence  and  this  question,  to  the  author's  knowledge, 
has  no  satisfactory  answer.  When  decorrelation  implies  statistical  indepen¬ 
dence,  there  is  no  specular  part  of  the  incoherent  power  because  the  second 
moment  of  the  scattered  field  behaves  asymptotically  (as  |  Ar ^ |  *•  “°)  exactly  like 
the  square  of  the  mean  scattered  field.  When  decorrelation  does  not  imply 
statistical  independence,  this  is  not  the  case  because  no  two  poinLs  on  the 
surface,  even  as  lAr^l-*00  ,  scatter  statistically  independent  of  each  oilier. 
That  is,  the  second  moment  of  the  scattered  field  does  not  approach  the  square 
of  the  mean  field  even  as  l^r^.l'*’00  .  As  a  final  word  of  caution,  one  should 
not  confuse  functional  dependence  with  statistical  dependence;  they  are,  in 
general,  completely  unrelated  mathematical  concepts. 

The  appearance  of  the  specular  term  in  (IS)  coupled  with  the  lack  of  any 
experimental  scattering  data  showing  this  effect  (to  the  author’s  limited 
knowledge)  suggest  that  decorrelation  does  indeed  imply  statistical  indepen¬ 
dence  for  real  surfaces.  However,  the  purpose  of  the  present  analysis  is  not 
to  address  this  question  but  rather  to  consider  the  theoretical  scattering 
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implications  of  statistically  dependant  surfaces.  It  is  therefore  beneficial 
to  illustrate  the  strength  of  the  specular  part  of  the  incoherent  scattered 
power  using  a  relatively  simple  example.  Such  is  the  intent  of  the  next 
section. 

Nume  r i ca 1  Examp 1 e 

In  order  to  keep  this  example  as  simple  as  possible,  the  surface  is 
stipulated  to  be  very  gently  undulating.  More  specifically,  the  surface  height, 
spectrum  Is  postulated  to  contain  no  spatial  frequencies  the  order  of  or 
larger  than  the  electromagnetic  wavenumber  kQ  and  the  variances  of  all  orders 
of  surface  height  derivatives  are  taken  to  be  arbitrarily  small.  Under  these 
conditions  the  physical  optics  approximation  for  the  current  induced  on  the 
surface  Is  valid,  l.e. 

5sC?)  -  2fi(r)  x  H1 (r)  (16) 

where  n(r)  ■  (-(3£/3x)x  -  0£/3y)y  +  2]/ [1  +  (3i;/3x)2  +  OcP y)2]  ^  is  the 
upward  directed  unit  normal  to  the  surface.  Substituting  (16)  in  (2),  taking 
the  stochastic  Fourier  transform  of  this  result,  and  then  separating  the  r  - 
dependence  off  as  in  (12)  yields 

jq(k,K)  -  2Ho(2x)i>”6(k-kiz)[c^(K)  -jC^'(Klx)6(K2y)  -  jcV  (kly)6(K,x)] 

(17) 

where  S„“l  +  and  the  K's  are  transform  variables  associated  with  the 

00  OO 

following  stochastic  variables  (with  p  ■  x  or  y); 

K  +-*■  V£,V2t,"  • 

K,  ♦+  3;/3p 

lp 

♦*  3t/3P  ,v2c  ,  V3C  ,  ••• 

The  constants  C.*1  -  q*(pxfi)  with  both  p  and  q  ranging  independently  over  x,  y, 

P 

and  z  are  related  to  the  polarization  of  the  Incident  field.  Substituting 
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07)  in  the  second  or  specular  term  In  (15),  denoted  as  Var(ty)  ,  yielJs  the 
following  result: 


Var(E8)  «  jr|"(2Tt)i  A6(ka  -k.  )  (2H  !)  f(l-k  -i)  (C  V  +  (l-k  -y)(i.y) 
“8  st  It  0  L  8  Z  S  z 

•<p,(k  -k  ,-k  +k .  ;(Ar  (■“)  -  p  (k  -k,  )p  (-k  +k .  )) 
^  2  sz  iz*  sz  iz’1  t1  rl  sz  iz  1  sz  iz  j 


0  ~ 

r 


08) 


where  and  are  now  the  Joint  two  point  and  marginal  c 


Imi  jctcrist  k  iunc- 


tlons,  respectively,  of  the  stochastic  height  only.  Normalizing  (18)  by  the 

i  1 2 

squared  magnitude  of  the  field  scattered  by  a  flat  surface,  |£°[  ,  and  noting 

that  specularity  implies  k  ■  -k  leads  to  the  following; 

sz  lz 


Var (Ea) 

_ _  s 

i*;!2 


P(-2kiz*2kiz;l^t,“*)  -P(-2kiz)P(2kiz) 


(19) 


The  squared  mean  scattered  field,  similarly  normalized,  is  given  by 


p(-2klz)p(2kiz) 


(20) 


The  specular  part  of  the  incoherent  power,  as  given  in  (19),  should  be 
compared  to  the  first  term  in  (15)  in  order  to  determine  which  is  dominant. 
However,  the  first  term  in  (15)  will  depend,  in  the  physical  optics  approxi¬ 
mation,  primarily  upon  the  slope  statistics.  Since  (19)  is  governed  by  the 
height  statistics,  it  is  possible  for  these  two  terms  to  exist  in  almost  any 
ratio.  A  more  meaningful  assessment  of  (19)  results  from  comparing  it  to 
(20).  Such  a  comparison,  in  essence,  relates  the  specular  portion  of  the  in¬ 
coherent  power  to  the  coherent  power  scattered  from  the  surface. 

For  this  example  the  Joint  probability  density  function  is  taken  to  be 
exponential  (Lennon  and  Papa,  1981],  e.g. 
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p(c,c' ;  I  Art  |-<*>)  - 


27l<^2> 


exp 


/?>/3 


(2)  ) 


2 

where  <£  >  is  the  height  variance.  The  joint  characteristic  fund  itm  can , 
with  the  aid  of  Fourier  transform  tables  [Campbell  and  Foster,  1961],  be 
shown  to  have  the  following  form; 


p(k,k' ;  | Ar 1 1 ■<*>) 


2  2 

(k  +k'  ) 


•3/2 


(22) 


The  marginal  characteristic  function  may  be  obtained  from  (22)  by  taking  k'=0. 
Substituting  these  results  in  (19)  and  (20)  yields 


Var(Es) 

mi2 


o  -3/2  ?  -3 

[1  +  8r/3]  -  [1  +  4£  / 3] 


(23) 


and 


\<*a>\2  2  "3 

- ~--[l+4Z^/3]  (24) 

ie;i 


where  E«*ko^C^>  coa0  is  the  Rayleigh  parameter  and  0  is  the  angle  oi  inci¬ 
dence  measured  from  the  z-axis.  Figure  1  illustrates  Lite  variation  of 
Var(E_)  /  | E°  J 2  (the  incoherent  specular  power)  and  | <E  > 1 2/| E° | “  (ihe  coher- 

“  fl  /  fl  b  *  S 

ent  power)  as  a  function  of  the  Rayleigh  parameter. 

When  the  Rayleigh  parameter  is  small,  the  incoherent  specular  power  is 
considerably  smaller  than  the  coherent  power  aa  it  should  be.  In  the  neigh¬ 
borhood  of  Z ■ 1  the  Incoherent  specular  power  is  equal  to  the  coherent  power 

and  for  E>1  It  greatly  exceeds  the  coherent  power.  For  large  Rayleigh  pararn- 

_3 

eter,  the  incoherent  specular  power  behaves  asymptotically  like  E  while  the 
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coherent  power  is  proportional  to  E”^. 

It  is  obvious  from  the  results  in  Figure  1  that  the  incoherent  specular 
power  cannot  be  ignored.  Since  the  total  power  scattered  from  the  surface 
must  be  in  agreement  with  the  incident  power,  the  existence  of  the  coherent 
specular  power  must  be  at  the  expense  of  the  diffuse  part  of  the  incoherent 
power,  i.e.  the  first  term  in  (15).  That  is,  the  Rayleigh  parameter  may  have 
to  be  relatively  large  before  the  diffuse  part  of  the  incoherent  power  is  sig¬ 
nificant  relative  to  the  incoherent  specular  and  the  coherent  powers.  Thus, 
when  comparing  the  diffuse  power  scattered  from  statistically  independent  and 
dependent  surfaces  (as  |Art|'*’fi0)  with  comparable  height  variances,  one  should 
expect  to  see  a  significant  difference  for  small  to  moderate  values  of  the 
Rayleigh  parameter.  For  very  large  values  of  E  this  difference  will  disap¬ 
pear  and  the  asymptotic  theories  of  Barrlck  [1968]  or  Sancer  [1969)  will  be 
valid.  For  statistically  dependent  surfaces.  Figure  1  illustrates  that  the 
validity  criterion  for  these  asymptotic  theories  should  be  based  upon  a  vanish 
lngly  small  value  of  incoherent  specular  power  rather  than  a  negligible  coher¬ 
ent  power. 

The  reader  is  reminded  that  the  purpose  of  this  paper  is  to  analytically 
describe  the  scattering  from  a  surface  for  which  decorrelation  does  not  imply 
statistical  independence.  This  paper  does  not  advocate  the  existence  of  such 
surfaces  in  the  real  world;  only  measured  statistics  can  answer  this  question. 
However,  these  results  do  show  how  electromagnetic  scattering  measurements  can 
be  used  to  identify  such  surfaces  provided  that  great  care  is  taken  in  separat 
lng  the  coherent  from  the  incoherent  scattered  power.  Finally,  these  analy¬ 
tical  results  also  show  that  extreme  caution  should  be  exercised  in  fitting 
statistically  dependent  functional  forms  for  the  probability  densities  to 
limited  histogram  measurements  if  one  is  interested  in  conditions  correspond¬ 
ing  to  1  $  £  S  10  . 
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Figure  Captions 

Figure  1.  A  comparison  of  the  coherent  and  incoherent  specular  powers  scat¬ 


tered  by  a  perfectly  conducting,  exponentially  distributed,  gently  undu¬ 
lating  ,  random  surface.  Both  powers  are  normalized  by  a  common  factor. 
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Abstract 

This  paper  derives  rigorous  results  pertaining  to  the  validity  of  the 
far  field  approximation  for  scattering  from  randomly  rough,  perfectly  conduct¬ 
ing  surfaces  having  arbitrary  statistics.  The  methodology  employs  the  stochas¬ 
tic  Fourier  transform  of  the  current  induced  on  the  infinite  surface  by  either 
a  bounded  or  unbounded  incident  plane  wave.  The  results  are  general  in  that 
no  approximate  simplifying  forms  for  the  current  are  employed.  Exact  expres¬ 
sions  are  obtained  for  the  mean  and  variance  of  the  scattered  field  for  unbounded 
illumination  and  they  are  compared  to  the  far  field  approximations  to  illustrate 
how  the  latter  simplifications  fail  in  this  limit.  Some  of  the  pitfalls  of  the 
far  field  approximation  in  the  case  of  beam  illumination  are  discussed.  When 
the  incident  plane  wave  is  bounded,  the  conventional  far  field  form  for  the 
mean  scattered  field  can  be  rigorously  derived  for  arbitrary  surfaces  pro¬ 
vided  the  cross  sectional  area  of  the  incident  beam  is  large  compared  to  the 
square  of  the  electromagnetic  wavelength.  The  conventional  far  field  result 
for  the  variance  of  the  scattered  field  is  shown  to  require  the  additional 
stipulation  that  the  cross  sectional  atea  of  the  incident  beam  contains  many 
decorrelation  intervals  of  the  surface  roughness.  The  results  obtained  herein 
are  important  because  they  hold  for  arbitrary  surface  statistics.  Whereas 
they  appear  to  duplicate  previous  results,  it  must  be  remembered  that  the 
earlier  results  were  only  valid  for  a  special  class  of  surface  statistics,  i.e. 
surfaces  for  which  single  scattering  theory  holds. 
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Introduction  and  Summary 

One  of  the  most  useful  descriptions  of  a  deterministic  body  scattering 

(or  radiating)  electromagnetic  energy  is  the  so-called  far  field  approximation. 

The  general  formula  for  the  electric  field  scattered  by  a  perfectly  conduct- 

-¥~ 

ing  body  and  evaluated  at  a  point  Rq  exterior  to  the  body  is  as  follows; 


V  x  V 
o  o 


/ 


J  (r)G( I R  -r I )  dS 
s  o 


(1) 


where  r)  =  /y  7e  is  the  characteristic  impedance  of  free  space,  k  =  2 tt / X 
o  o  o  o  o 

is  the  electromagnetic  wavenumber,  G  is  free  space  Green’s  function 

■  -*■  *  >  -fi 

exp(  -  j  kQ  |Ro  -  r  |  )/4it|Ro  -  r  I  ,  and  Jg  is  the  current  density  on  the  surface 
S  of  the  body.  If  the  point  Rq  is  sufficiently  far  removed  from  the  body 
that  (2)  may  be  approximated  by  the  following  expression  [Van  Bladel,  1964] 


E  (R  )  s 
s  o 


j  k  n  G(R  )  k 
o  o  os 


X  k 


‘fa 


)  exp  (j  kokg« 


r)  dS 


(2) 


A  -t- 

where  kg  Is  a  unit  vector  specifying  the  direction  of  the  point  Rq  ,  then 
(2)  is  essentially  the  definition  of  the  far  field  scattered  by  the  body.  The 
results  of  extensive  computations  by  Hansen  and  Bailin  [1959]  have  provided 
much  insight  into  the  magnitude  of  Rq  necessary  to  accurately  replace  (1) 
by  (2)  for  deterministic  currents.  Consequently,  the  range  of  validity  of  the 
far  field  approximation  is  reasonably  well  understood  for  deterministic  prob¬ 
lems. 

When  the  geometry  of  the  scattering  body  varies  on  a  sample-by-sample 
basis  so  that  it  i3  more  conveniently  described  in  a  probabilistic  manner, 
the  concept  and  implications  of  a  far  field  are  less  well  understood.  The  lack 
of  clarity  for  this  situation  stems  from  two  difficulties  which  are  not 
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encountered  in  the  deterministic  case.  First,  for  the  deterministic  geometry 
one  may  postulate  certain  cannonical  current  distributions,  compute  Eg  from 
(1)  and  (2),  and  obtain  very  useful  information  on  the  range  of  validity  of 
the  far  field  approximation.  However,  for  the  random  geometry  this  cannot  be 
done  because  the  dependency  of  the  current  upon  the  random  geometry  is  criti¬ 
cal  and  it  must  be  obtained  (not  postulated)  from  the  boundary  equations. 

The  second  difficulty  stems  from  the  need,  for  the  random  geometry  problem,  to 
know  not  only  Eg  but  all  or  a  portion  of  its  statistical  moments  and  these 
are  not  as  simply  related  as  the  scattered  field  and  power  in  the  determinis¬ 
tic  case. 

A  few  studies  have  obtained  useful  results  on  the  random  geometry-far 
field  problem  [Barrick,  1965],  [Miller,  1982],  [Fung  and  Eom,  1982];  however, 
they  are  limited  by  the  use  of  an  approximate  form  for  the  current  and  only 
the  Inclusion  of  a  quadratic  phase  term  in  (2).  What  is  really  needed  is  some 
general  insight  into  the  far-field  approximation,  as  it  applies  to  the  random 
geometry  situation,  which  is  not  limited  by  restrictive  simplifications.  Such 
is  the  goal  of  this  paper. 

The  first  problem  to  be  considered  is  that  of  a  randomly  rough,  perfectly 
conducting  surface  which  is  infinite  in  horizontal  extent  and  illuminated  by 
an  unbounded  plane  wave.  This  infinite  geometry/ illumination  case  is  selected 
as  the  starting  point  because  the  current  on  the  surface  assumes  a  somewhat 
simplified  form  and  this  permits  the  derivation  of  some  useful  exact  results. 
Although  the  far  field  assumption  is  not  applicable  to  this  problem,  results 
obtained  from  this  approximation  are  reviewed  in  order  to  illustrate  the  func¬ 
tional  dependencies  predicted  by  it.  The  mean  and  the  variance  of  scattered 
field  are  next  obtained  using  the  exact  expression  for  the  scattered  field, 
i.e.  (1).  The  mean  scattered  field  is  shown  to  be  a  plane  wave  propagating 
away  from  the  surface  in  the  specular  direction.  The  amplitude  of  this  plane 
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wave  is  reduced  relative  to  the  amplitude  of  the  incident  plane  wave  by  the 
effects  of  the  roughness  on  the  surface.  The  polarization  of  the  scattered 
plane  wave  is  dependent  upon  the  multiple  scattering  processes  occurring  on 
the  surface  and  this  shows  that  the  mean  scattered  field,  contrary  to  previous 
statements  [Moore  and  Williams,  1957],  is  governed  by  the  laws  of  diffraction 
and  not  reflection. 

The  exact  expression  for  the  variance  of  the  scattered  field  is  found  to 
be  considerably  more  complicated  than  the  mean  scattered  field  relationship. 
However,  the  exact  expression  clearly  shows  that  the  far  field  approximation 
is  not  applicable  when  the  surface  and  the  illumination  are  unbounded  even 
though  the  point  of  observation  is  taken  to  infinity. 

The  results  obtained  for  the  unbounded  geometry /illumination  problem  are 
important  from  a  theoretical  point  of  view  but  they  have  little  significance 
relative  to  the  practical  situation  wherein  the  illumination  is  bounded,  i.e. 
a  beam  plane  wave.  When  the  incident  illumination  is  bounded,  the  problem 

becomes  difficult  because  the  surface  current  is  more  complicated.  However, 

2 

it  is  shown  that  when  the  illuminated  area  is  very  large  compared  to  , 
the  current  appropriate  to  the  unbounded  illumination  case  in  conjunction 
with  the  bound  on  the  incident  illumination  may  be  used  with  little  error. 

This  result  for  the  current  is  then  used  in  (1)  for  the  purpose  of  com¬ 
puting  the  mean  and  variance  of  the  scattered  field.  When  multiple  scattering 
on  the  surface  is  negligible,  this  procedure  produces  an  exact  result.  When 
multiple  scattering  is  important,  the  procedure  leads  to  an  underestimation  of 
the  true  illuminated  area  on  the  surface,  i.e.  the  true  illuminated  area  is 

the  area  over  which  the  current  is  nonzero.  However,  as  long  as  the  geometric 

2 

cross  sectional  area  of  the  incident  beam  is  large  compared  to  ,  the  error 
in  this  procedure  is  small.  Proceeding  with  this  approach,  i-  is  possible  to 
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show  in  a  very  clear  manner  why  the  far  field  approximation  for  the  scattered 
field  and  the  case  of  a  randomly  rough  scattering  surface  is  not  valid  as  the 
illuminated  area  becomes  infinite.  Finally,  conditions  are  derived  under 
which  the  far  field  simplification  is  valid  for  arbitrary  surfaces. 

While  the  far  field  approximation  has  been  used  extensively  in  randomly 
rough  surface  scattering  problems,  justification  of  the  approximation  has  here¬ 
tofore  been  limited  at  best.  The  results  presented  in  this  paper  provide  a 
general  foundation  for  the  approximation  and  yield  insight  into  why  it  fails 
for  an  infinite  surface. 

Background 

The  surface  described  by  z  =  £(x,y)  is  assumed  to  be  infinite  in  the 
x  and  y  directions  and  for  z  >  £(x, y)  the  medium  is  taken  to  be  free  space 
while  for  z  <  £(x,y)  the  medium  is  perfectly  conducting.  The  stochastic 
surface  height  £(x,y)  comprises  a  zero  mean,  statistically  homogeneous  process 
with  the  mean  surface  corresponding  to  the  •  =  0  plane.  A  plane  wave  illumi¬ 
nates  the  surface  and  its  spatial  support  in  the  x-y  plane  will  be  taken  to  be 
either  bounded,  i.e.  a  beam  plane  wave,  or  infinite.  The  explicit  form  of  the 
incident  plane  wave  is  as  follows; 

E1  =  Eq  exp(-  j  £t*r)  9  (rt>  (3) 

where  k^  =  kQk^  =  k^  +  k^  z  and  k^  specifies  its  direction  of  propagation. 

c  z 

The  function  0(rt>  represents  the  support  of  the  incident  plane  and  its 
dependence  upon  r^  *  xx  +  yy  signifies  that  it  is  independent  of  the  z- 
dimension. 

/v 

The  mean  field  scattered  in  the  direction  kg  and  having  a  polarization 
e  is,  in  the  far  field  approximation,  given  by  the  following; 


<E  *e>  *  T  e  *  k  x  k  x  1  (3  (r)  exp  (j  k  £  +  j  k  *r  )  db)> 
s  s  s  I's  s  s  t 

J  2  L 


S 


S 


where  k  +k  z  =  k  k  and  T  «  i  k  n  G(R  )  .  To  simplify  matters  somewhat, 
ss  os  Joo  o 

t  z 

the  integration  over  the  rough  surface  is  transformed  into  an  integration  over 

the  z  =  0  plane  through  dS  =  ^1  +  (3£/3x)“  +  (3£/3y)^  drt  where  drt  =  dxdy 

and  3£/3x  and  3C/3y  are  the  surface  slopes  in  the  x  and  y-directions.  Thus, 

*>■  •> 

defining  the  modified  current  J(r)  as 


J(r) 


ViT 


OC/3x)2  +  (3C/3y)2  Js(r)  , 


the  far  field  approximation  for  <6^  •  e>  becomes 


:E  •  e>  =  T  e*k  *  k  x 
s  s  s 


/O(r)  exp  (j  kg  £))>exp(jks  *rt)  drt 
z  t  L 


(4) 


The  average  in  (4)  may  be  explicitly  written  as  follows; 


<J(r)exp(jk  C) 


>  =  j ‘  \f  J^pl^^,V2,,-.)exp(jks  C)dCdVC---  (5) 


where  is  symbolic  for  the  surface  slopes,  7  £  is  symbolic  for  the  sur¬ 
face  curvatures,  etc.  The  function  in  (5)  is  the  single  point  joint 

probability  density  function  for  C»  V£,  etc.  The  integrations  in  (5) 

may  also  be  written  as  convolutions  of  transforms  with  respect  to  the  stochas¬ 
tic  variables,  i.e. 


<•> 


(2  tt> 


d^t’S.’k2’*  *  *  >  ^l^s  "k1»_k2’"  ’  ■*  dkld^2 


(6) 


where  s_  *  l  i  and 

CO 

i-1 
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J=j'*j  J(r)  exp  ^jk1C  +  j  Vnr  .  kn+1^  dC  dVf,-  • 

**  n=l 

1  =/“/  eXp(jkl?  +  j  2  ^  ‘Vl) 


(7a) 


dCd7C--'  (7b) 


With  K.  signifying  the  dependence  upon  ^2 ’^3’"*  ’  a  su^st- i tution  of  (6)  into 
(4)  yields  the  following; 


<E  «e>  =  T  e  •  k  xk  x  - i— 

s  s  s  s 

(2tt )  ° 


HI  J (rt  ,k^  ,K)  (k^  -k^,-K 


K)exp(jks  'r^dr^dk^dk 
z  t 

(8) 


It  should  be  noted  that  is  the  single  point  joint  characteristic  function 

2 

for  the  random  variables  £,V£,V  £,**'  . 

In  order  to  proceed  further  one  must  have  some  knowledge  of  the  stochastic 

Fourier  transform  of  the  modified  current,  i.e.  Jfr^.k^.K)  .  In  particular, 

■>  -> 

the  variation  of  J  with  r  is  essential  to  determining  the  directional 

characteristics  of  the  mean  scattered  field,  i.e.  the  dependence  of  <Eg  •  e> 

upon  k  .  The  only  case  for  which  the  exact  r  -dependence  cf  J  is  known 
St  C 

is  when  the  incident  plane  wave  is  unbounded  or  0(rt>  =  1  for  all  r  .  When 

this  occurs,  it  has  been  shown  [Brown,  1982b]  that  J  assumes  the  following 

form; 

J(rt»ki»K)  =  j(klsK)  exp  (-jki  *rt)  (9) 

Substituting  (9)  in  (8)  and  taking  e  orthogonal  to  k  in  the  specular 

s 

direction  yields 

<Eg*e>  -  (2Ti)2r<5(kg  -1^  )  I  I  e  •  jCk^K)  4>1  (kg  -k^-K)  dkjdK/^ir)  °°  (10) 

t  t  J  J  z 


This  result  indicates  that  the  mean  scattered  field  is  a  spherical  wave  modified 
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by  a  product  of  two  amplitude  weighting  factors  -  the  6-function  and  the  inte¬ 
gral  term.  The  integral  term  is  entirely  a  consequence  of  the  surface  rough- 


ness.  The  6(k  -k  )  factor,  loosely  interpreted,  shows  that  the  scattered 
St  it 


A  A 


field  is  nonzero  only  in  the  direction  k  =k  .  Furthermore,  since  the  point 

st  it 

of  observation  is  above  the  surface  the  condition  k  51  -k.  must  also  be 

sz  iz 

satisfied.  The  form  of  (10)  is  certainly  less  than  satisfying  and,  in  fact, 
only  has  meaning  when  compared  to  the  field  "scattered"  by  a  flat  plane  in  the 
far  field  approximation  [Brown,  1982a].  In  short,  (10)  cannot  be  taken  liter¬ 
ally  but,  instead,  must  be  interpreted. 

The  second  moment  of  Eg  • e  ,  in  the  far  field  approximation,  is  obtained 
by  forming  the  product  of  (2)  and  its  complex  conjugate  and  then  averaging 
over  all  stochastic  variables.  The  averaging  operation  requires  use  of  the 
two  point  joint  probability  density  function,  i.e.  p  (C ,£_,V!; .,VC  ,• • • ;r  -r  ), 

^  1  fc  J-  t 

which  depends  not  only  upon  the  random  surface  characteristics  at  r  and 


r  but  also  on  the  vector  distance  r  -r 

2 


-r  .  Since  E  is  polarized  orthog- 

H  c2  s 

onal  to  k  in  the  far  field  approximation,  e  is  also  chosen  to  be  orthogonal 


to  kg  and  the  second  moment,  using  the  far  field  approximation,  can  be  written 
as  follows; 

;|E,‘*|2>  -  |r]2Jj[l(?1)'SHJ*(r2)'SJp2(S1,i;2,VCJ,VC2,---;rt  -rt  )exp[jks* (r^-r^ ] 


Converting  the  Integrations  over  the  stochastic  variables  to  convolutions  in 
transform  space  and  using  (9) ,  the  above  result  can  be  manipulated  into  the 
following  form; 
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<  Ie  *e  I  2>  *  lim  *"•. —  j  • 

'  8  '  A-  (2»)2s*  J 


Ha** 


K^’eltTVk^-K^-S] 


*2«s  -kj^.-ks  -k2»-K1,-K2;Art)exp^J  (kg  )*ArtJ  dArtdk1dk2dK1dK2 


where  A  is  Che  illuminated  area  and 


<£2(ki»k2»^i’^2;A^t)  =  ^****|'p2(?i»52’3?l’3?2;A^t^exp[^kir,l  +  k2Zjl 

+  ^•3^  +  K2*3?2)J  d;id?2d(3c1)d(3ra2)  (12) 

is  the  two  point  joint  characteristic  function  with  3£^  and  3^2  symbolic 
for  all  orders  of  surface  height  derivatives  Vn£,n-1,2 ,  *  * ‘  .  Assuming  that 
the  surface  statistics  are  such  that  decorrelation  implies  statistical  indepen¬ 


dence  yields 

_  lim  P2(C1,C2,3C1,3^2;A:?t)  -  0 

lArJ-x” 

This  condition  implies  that  <J?2  will  have  infinite  support  with  respect  to 
Art  and  so  the  transform  of  cj>.,  contained  in  (11)  will  yield  a  6-function 
dependence  upon  kg  -k^  .  This  problem  can  be  overcome  by  using  a  function 
in  the  integrand  of  (11)  which  has  a  finite  support  with  respect  to  Art  . 
Thus,  by  subtracting  and  adding  ct^O,*;00)  to  the  integrand  of  (11)  yields 


<|E  *e 1 2>  -  lim  ■  A^F-1 — 


|  E  -e |  >  «  lim 
8  A-«° 


|p  J’  •  'J l  j  (kx .  Kj_ )  *e  1 [  j ’ *  (-k2  ,-£,)•£]  j< 


^2^8  ~kl’  ks  k2’-Ki»_K2;Art 
z  z 


-  <3>2 (ks  -k1,-kg  -k2,-Klt-K2;ooA  expjj  (kg  -k^  ) •ArldAr^^d^d^c 
sz  sz  J  L  t  t  J 

+  t(VVj  "/ ^(kl-Sl)-ai(T*(-k2--iE2)-81 


*2<k8  -klf-ks  -kj.-l^.-^t-Jdkjdk^d^ 
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Because  of  (13),  may  be  rewritten  as  a  product  of  marginal 

characteristic  functions.  Furthermore,  since  the  surface  roughness  is  a  zero 
mean  process,  the  marginal  characteristic  functions  will  be  real  and  symmetric. 

Returning  to  the  integral  definition  of  6(k  -k  )  ,  it  may  be  shown  that 

2  8t  t 

1 6 (k  -k  )|  should  be  Interpreted  as  follows; 

St  *t 


|«£.  -*i  )\2  ■ lim  A  6(ks  ) 

1  St  Xt  1  A-m»  ®t 


Thus,  the  last  term  in  (14)  is  recognized  to  be  the  magnitude  squared  of  the 
mean  field,  i.e. 


|<ts-a>|2  -  a.  -4^,  set.  >1  f[e-l<h’ii>*i<k,  -V-i<2>dkidlti 
*“  (21.)  *  c  c  IJJ  * 


and  so  the  variance  of  the  scattered  field  is  as  follows; 


(15) 


Var (E  »e)  -  lim 
3  A-*» 


(2ir) 


K2)*e) 


*2(k8  -klt-k  -k2»*^i‘_K2;ks  >dk1dk2dK1dK2  (16) 


z  z 

1 2  _  /_  ,2  z 


t  t 


where  |r|  -  »  <*>2(Art)  *  ^(Ar^  -  <X>2(»)  ,  and 

•2<--»«s  *  “  J  |  ’*  “<**2^*  ’*  *“)}exP[j<ks  -ki  (16a> 


The  Important  point  to  note  about  (16)  is  that  *-he  directional  characteristics 
of  Var($s*£)  are  determined  essentially  by  th-  r  transform  of  «t>2 

A  A 

evaluated  in  the  direction  k  -k  .  That  is,  the  far  field  approximation 

8t  1t 

predicts  that  the  variance  of  the  scattered  field  is  proportional  to  only  one 
Fourier  consonant  of  <j>2  • 
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The  results  obtained  in  this  section  are  not  intended  to  justify  the  uBe 
of  the  far  field  approximation  when  the  surface  extent  and  the  illumination 
are  both  unbounded.  Rather,  this  section  is  intended  to  serve  two  other  pur¬ 
poses.  First,  the  technique  of  treating  averages  as  convolutions  in  stochas¬ 
tic  transform  space  has  been  introduced  and  used  for  the  purpose  of  familiari¬ 
zation.  Second  and  more  important  is  the  fact  that  the  only  approximation  con¬ 
tained  in  this  section's  results  is  the  far  field  approximation.  Thus,  it  is 
now  possible  to  obtain  similar  results  from  the  exact  expression  for  the  scat¬ 
tered  field  and  compare  them  with  (IS)  and  (16)  without  having  the  comparison 
tainted  by  the  use  of  an  approximate  current.  That  is,  any  differences  between 
(15)  and  (16)  and  corresponding  exact  results  will  be  solely  a  consequence  of 
the  far  field  approximation. 

Exact  Results  for  Unbounded  Illumination 

The  exact  mean  scattered  field  can  be  obtained  without  much  difficulty 

from  (1);  however,  this  is  not  necessary  since  it  can  also  be  derived  from 

<|?  *e|^>  ,  apart  from  a  sign, by  the  same  approach  as  demonstrated  above.  The 
8 

^  2 

derivation  of  a  simplified  expression  for  <|Ea*ep>  from  (1)  is  also  not 
difficult  but  it  is  involved  and  normally  it  would  be  better  to  relegate  it 
to  an  appendix.  However,  it  does  contain  a  very  important  simplifying  step 
and  so  it  will  be  included  in  the  main  text.  To  avoid  a  minor  difficulty 
[Brown,  1982a],  surfaces  for  which  decorrelation  implies  statistical  inde¬ 
pendence  will  be  assumed. 

Using  (1)  and  the  technique  of  converting  averages  into  convolutions  in 
transform  space  yleldB  the  following; 
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<(Es*e|2>  «  [n0/k0]2j%’“J'  ^  ^2  ^2  (-k1»“k2»_ki>-k2;^t^t 


2s. 


•  drt  dk1dk2dK1dK2  /<2ir) 


(17) 


where  Ar  »  r  -r  and 
t  t2 


F{  }1  -  (27T)"1  e  •  Vo  x  VQ  x 


J  J  J(^,k1-B1,K1)G(Ro  -^fZo-C1)exp(JB151)dC1d81 

F{  )2  *  (2ir)_1  e  •  Vq  x  VQ  x  J  J  J*(r^  .k^  ,K2)G*(Rq  -r^  ,2q  -  C2)exp(jB2C2)dC2d£ 


and  R  «  X  x  +  Y  ?  is  the  horizontal  distance  from  the  origin  on  the  surface  to 
o^  o  o 

the  point  of  observation.  The  tilde  symbol  over  J  and  J*  denotes  the 
Fourier  transform  with  respect  to  the  stochastic  variables.  Since  both  the 
surface  and  the  illumination  are  unbounded,  (9)  may  be  used  along  with  the 
following; 

J*(rt  ,k2-B2,K2)  -  {  J(^,-k2+e2,-K2)}* 


Thus,  (17)  reduces  to  the  following  form; 


<  |  E  -S 
1  s 


1 2>-<W2/--/  {8-,o^oX^<kl-ei-Jl><!(*o*'t1,VCl)expl-^lt''t1'fj8lCl1} 


•  {«-W^<-k2^2--k2)C*<kot-?t2'V42,“['(+:|i‘it'?t2  +162C2>} 


<r>2(-k1,-k2,-K1,-it2;A^t)drt^dr^dC1d?2dB1dB2dk1dk2dK1dK2/(2Ti) 


2  (soo+1) 


(18) 
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The  r  -integration  involves  only  the  functions  C  and  <J>,  .  Using  the  shift 
C1  1 
theorem  for  Fourier  transforms  and  the  two-dimensional  Fourier  transform  of 


G  ,  the  r^-integration  may  be  expressed  as  follows; 


where  the  tilde  over  <t>2  denotes  the  transform  from  Art  to  kt  and 

(k  -k  )2  ■  (k  -it  )*(k  -it  )  .  Since  the  point  of  observation  (X  ,Y  ,Z  ) 
t  if  t  it  t  o  o  o 

is  to  be  above  any  realization  of  the  rough  surface,  Iz  -£. I  *  Z  -L,  .  Sub- 

o  JL  o  1 

stituting  these  results  in  (18)  and  performing  the  Cj- integration  gives  rise 
to  a  delta  function,  i.e.  6^(5^  +  ^k2  -  (k^  -k^2*  ^  ,  which  makes  the 
integration  trivial.  Thus,  remembering  that  Zq  must  be  above  the  surface, 
(18)  reduces  to  the  following; 


<n  /  k  ) 

O  O 


e*VV^kl+K*Kl)  expt-jtcz0+jkt- 


~r  )]/(-j2ie) 

t  *2 


•  “»w£i  •*„  )}{s-7ox7o4*<-k,«2.-£2)o*(?<i  -;t  .zo-i2> 


exp[+jk. 


ft2+jB2C2] 


;-kt)d?2dB2d^t  dktdk1dk2dK1dK2  f  (2tt) 


2sm+3 


2  2. 
where  ic  “  [k  -  (k.  -k  )  ] 
o  it  t 

dimensional  Fourier  transform  of 


The  r  -integration  in  (20)  yields  the  two- 
C2 

G*  which  contains  an  exponential  term  whose 


argument  is  proportional  to  |7.o-52|  *  Dropping  the  magnitude  sign  since  Zq 
is  to  always  be  greater  than  ?2  and  performing  the  ^“integration  in  (20) 
yields  a  6-function  which  simplifies  the  B2-integration.  The  result  is  as 
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follows; 


<iV«i2>  -  <w2// [ Mi  '/(-i2E)} 

’  {S*  W  j(k-k2+k,-K2)  exp [-jk  Zo-jkt-Ro  ]  j (-j2k) j* 


28^+2 


•  <t>2  (-Viv-k2,-lLL,-K.2;kt-ki  ^dk^dK^/tfir)  (21) 

2  -+  2  1/2 

where  k  *  [kQ  -  kfc  ]  .  Regularizing  the  integrand  in  (21)  so  that  it  does 

not  contain  a  6(kt~ki  )  dependence  entails  subtracting  and  adding 
prior  to  transforming  from  Ar^  to  k^k^  •  Thus,  with 

r 

<t>2^'  ’*  ;kt~ki  *  "  I  ^2^* ’*  ;Art^  "<^>2^‘,‘;“^exp^^t_^i 
t  J  t 

and  the  fact  that  nay  be  written  as  a  product  of  marginal  charac¬ 

teristic  functions,  the  regularized  form  of  (21)  becomes 


<  t  Es  -S  |  > 


<W/////  {swi^'V*(Et'V){J'W»(¥'-|2»*|vv)‘ 


•  i2(-k1,-k2,-K1,-K2;kt-ki  )dktdkldk2dK1dK2/(2TT) 


2s  +2 


<vv2/ f J 


*  <*>!  (-k1+>/k q2-  k±2  ,  -K^  4^ (k2  -VkQ2-  k£2  ,  K2)dk1dk2dK1dK2  / (2n)  "  (22) 


^  ^  n  l/2 

where  Rq  -  Xo*  +  YQy  +  Zqz  ,  k  -  (k^  -k£)  ,  and 
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(22a) 


\|)(k  ,R )  ■  exp  [-  jk  Z 
to  o 


]  /(-j  2k) 


The  last  term  in  (22)  is  the  magnitude  squared  of  the  mean  scattered  field, 
|<E  *e>|  ,  and  the  following  operations  demonstrate  this  fact.  Since  (J>  is 
real  and  symmetric  in  all  transform  variables,  the  last  term  in  (22)  which  re¬ 
sults  from  regularizing  (21)  is  recognized  to  be  the  product  of  one  set  of 
integrals  and  its  complex  conjugate.  Thus,  it  may  be  expressed  as  follows; 


W  J  (k^wtfT  ,Ro)l<D1(-k1+Vk^k71 

t  '  it 


dkjdK- 


(23) 


The  function  ^(k^,^)  is  independent  of  k^  and  and  is  a  plane  wave.  The 
1/2  ■+  -> 

root  of  (k  -k  )  must  be  chosen  such  that  lp(k^  ,R  )  is  upward  traveling 
t  c  °  ^  ^ 

because  is  positive;  this  dictates  the  following  form  for  ip (k .  ,Rq)  ; 


^(k±  ,Ro)  ■  exp  T~  J  |ki  !  Zq  -  jk  -Ro  ] /(- j  2  !k±  |) 
t  L  z  t  t  ‘  z 


(24) 


where  the  magnitude  signs  on  k  are  used  because,  according  to  (3),  k.  <_  0  . 

Z  lz  " 

The  quantities  VQXV  x  j'J'  may  be  manipulated  using  vector  identities  into  the 

form 


VVJ*-kol>}+  w]  (25) 

By  expanding  the  vector  operations  in  (25)  with  the  aid  of  (24),  it  may  be 
shown  that 


7o  x  Vo  x  T  ij.  -  [k®  x  k®  x  l]ko  *8/(  ~  J2cos  ©i) 
where  9^  is  the  angle  of  incidence  measured  with  respect  to  the  z-axis. 
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kg  is  a  unit  vector  in  the  specular  direction,  i.e. 


k  =  <k,  +  lk,  I  z)/  +  lk,  I  z|  » 


it  '  *z 


it  .~iz 


and 


ip  =  exp  (  -  j  k  R  k  *k) 
rs  K  Joos  s 


Substituting  these  results  in  (23)  and  taking  e*kg  =  0  yields 


l<VH 


(26) 


and  it  can  be  shown  that  the  expression  inside  the  magnitude  brackets  is  iden- 

•f  A 

tical,  except  for  a  minus  sign,  to  the  average  of  Eg*e  obtained  from  (1). 
Substracting  (26)  from  (22)  gives  the  following  relationship  for  the  variance 
of  the  scattered  field; 


Var(Eg-e) 


'Wj  <ki+k»K1)^(kt,Ro) 


}{* 


•V  xV  xj(-k  +k,-K  ) 
o  o  z  z 


Mkt,R0 


>}' 


<t>2  (-kx  ,-k2  ,-K^  ,-K2  5kt-ki  )dkj.dk1dk2dK1dK2 


/(2tt) 


2  s  +2 
00 


(27) 


Unfortunately,  the  integrations  in  (26)  and  (27)  cannot  be  accomplished  in 
general  without  knowledge  of  jf,  <3^,  and  <t>2  . 

The  important  point  to  remember  about  (26)  and  (27)  is  that  they  are  both 
exact  results.  Furthermore,  it  is  not  necessary  for  the  point  of  observation 
Rq  to  be  infinitely  far  from  the  surface.  The  only  restriction  on  these 
results  is  that  the  point  of  observation  be  above  any  realization  of  the  surface. 


Comparison  of  Exact  and  Far  Field  Results 

Equations  (15)  and  (16)  are  results  obtained  using  only  the  far  field 
approximation  while  (26)  and  (27)  are  exact.  A  comparison  of  the  two  sets  will 
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show  how  che  far  field  approximation  breaks  down  when  both  Che  surface  and 
Che  illumination  are  unbounded. 

The  exact  result  for  the  mean  scattered  power,  i.e.  (26),  shows  that  the 
mean  scattered  field  is  a  plane  wave  propagating  in  the  specular  direction, 

A  A  A  A 

s  s 

k  »  k  and  k  «  -k.  .  Furthermore,  it  can  be  shown  that  (26)  yields 

s  s ,  S  1 

t  i  z  z 

the  correct  reflected  power  when  the  random  surface  goes  to  a  plane.  The  far 
field  approximation  for  <Es*e>  ,  i.e.  (10),  yields  a  spherical  wave  with  an 
amplitude  weighting  such  that  it  is  nonzero  only  in  the  specular  direction. 

In  the  limit  of  the  random  surface  approaching  a  plane,  the  far  field  approxi¬ 
mation  does  not  reduce  to  a  reflected  plane  wave  traveling  in  the  specular 
direction.  It  is  interesting  to  note  that  although  the  far  field  approximation 
does  not  produce  the  plane  wave  nature  of  the  scattered  field,  it  does  contain 
the  proper  dependence  upon  the  statistical  properties  of  the  surface  because 

in  the  specular  direction  k  «*  |k  |  . 

s  iz 

z 

Equation  (26)  does  not  provide  any  insight  into  the  polarization  of  the 

mean  scattered  field;  the  polarization  is  determined  entirely  by  the  transformed 

current  T(k^,lt^)  .  If  there  is  no  significant  multiple  scattering  occurring 

on  the  rough  surface  then  <Eg>  will  have  the  same  polarization  as  a  field 

reflected  from  a  perfectly  flat  surface.  In  this  case  <E^>  will  have  all  of 

the  characteristics  of  a  specularly  reflected  field  except  that  it  will  be 

attenuated  by  the  effects  of  the  surface  roughness.  If  multiple  scattering  is 

significant,  the  polarization  of  <Eg>  will  be  different  from  that  of  a  planar 

surface  reflected  field.  This  result  shows  that,  contrary  to  previous  thoughts 

[Moore  and  Williams,  1957],  the  mean  scattered  field  is  not  governed  by  the 

laws  of  reflection  if  multiple  scattering  on  the  surface  is  significant. 

Furthermore,  it  is  not  advisable  to  refer  to  the  mean  scattered  field  as  a 

specular  field  since  <Eg>  does  not  always  have  the  characteristics  of  a  true 
specularly  reflected  field. 
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According  to  (16)  the  variance  of  the  scattered  field  is  linearly  related 

to  the  illuminated  area  and  inversely  proportional  to  the  second  power  of  the 

distance  from  the  origin  on  the  surface  to  the  point  of  observation.  More 

importantly  however  is  the  fact  that  the  variance  is  determined  by  the  Fourier 

♦ 

transform  of  the  regularized  characteristic  function  evaluated  at  k  - k .  , 

st  4t 

i.e.  see  (16)  and  (16a).  This  is  a  direct  consequence  of  the  fact  that  if  one 
represents  the  scattered  field  in  terms  of  an  angular  spectrum  of  plane  waves 
[Clemmow,  1966],  the  far  field  approximation  leads  to  a  scattered  field  which 
is  determined  by  one  and  only  one  angle  or  direction  in  the  spectral  represen¬ 
tation.  In  terms  of  the  variance  of  the  scattered  field,  this  means  that  one 
obtains  a  result  which  is  determined  by  one  and  only  one  Fourier  component 
of  the  regularized  characteristic  function.  Contrast  this  result  with  the  exact 
result  as  given  by  (27)  in  which  the  variance  of  the  scattered  field  is  deter¬ 
mined  by  a  weighted  integral  of  all  Fourier  components  of  the  regularized 
characteristic  function.  In  fact,  the  very  presence  of  the  ^-integration  in 
(27)  shews  that  there  is  no  such  thing  as  a  far  field  for  the  case  of  unbounded 
illumination  regardless  of  how  far  the  point  of  observation  is  moved  away  from 
the  surface.  For  example,  by  moving  the  point  of  observation  away  from  the 
surface  one  merely  sees  less  of  the  surface  waves  propagating  along  the  surface. 

In  fact,  in  the  limit  of  Rq  -► 40  the  k^-integrat ion  in  (27)  is  essentially 

■+2  2  ■*  > 

limited  to  the  visible  range  (k  <k  )  because  of  the  iMk^,R  )  function; 

to  to 

but  Var(Eg*e)  is  still  determined  by  a  weighted  integral  of  all  Fourier 
components  of  over  the  visible  range  of  wavenumbers.  As  an  aside,  it 

should  be  noted  that  although  one  does  not  observe  any  surface  waves  as  R0_vo° 
they  still  influence  the  result  through  the  transformed  current  j  .  As  a 
final  point,  it  can  be  shown  that  (27)  goes  to  the  proper  limits  in  the  cases 
of  a  planar  surface  or  a  randomly  elevated  planar  surface  whereas  (16)  does  not. 
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It  would  be  very  beneficial  if  the  ^-integration  in  (27)  could  be  accom¬ 
plished  without  explicit  knowledge  of  either  J  or ^  •  However,  this  is  not 
possible  primarily  because  the  convergence  of  the  integral  is  highly  dependent 
upon  j  and  ^  •  This  can  be  seen  by  noting  that  if  the  j's  and  <t>2  arc  set  to 

■>  i  -  i  n 

unity,  the  k  -integration  is  exactly  equal  to  the  integral  of  »  [ e " V  *G(R  )  1 

t  o  o  o 

over  the  x-y  plane  and  this  integral  does  not  converge. 

It  is  certainly  not  surprising  that  the  far  field  approximation  fails  when 
the  incident  illumination  is  unbounded.  In  the  case  of  the  mean  scattered 
field,  it  has  been  shown  that  the  effects  of  the  surface  roughness  are  correctly 
predicted  by  the  far  field  approximation.  However,  for  the  variance  of  the 
scattered  field  there  is  a  very  fundamental  difference  between  the  far  field 
results  and  the  exact  results  and  so  the  far  field  approximation  is  essentially 
useless. 

Bounded  Illumination 

With  the  exception  of  certain  fundamental  theoretical  questions,  the  case 
of  bounded  incident  illumination  is  of  most  interest.  Remembering  that  the 
intent  here  is  to  obtain  general  results  which  do  not  depend  on  explicit  approxi¬ 
mate  forms  for  the  surface  current,  the  bounding  of  the  illumination  causes  a 
problem  because  the  separation  in  (9)  is  no  longer  valid.  In  order  to  under¬ 
stand  this  difficulty  and  investigate  appropriate  approximations,  it  is  neces¬ 
sary  to  turn  to  the  stochastic  Fourier  transformed  magnetic  field  integral 

— -W 

equation  for  the  modified  current  J  [Brown,  1982b].  This  equation  is  obtained 

CD 

by  multiplying  the  magnetic  field  integral  equation  by  exp[jk,t;  +  j  I  "k  . ,  *Vn^] 

1  n=l  n+i 

and  then  averaging  over  all  stochastic  surface  characteristics  Vn£ ,n=0 , 1 , •  •  •  . 

The  averages  are  converted  to  convolutions  in  stochastic  transform  space  and 
one  obtains  an  infinite  dimensional  integral  equation  of  the  first  kind  having 


che  following  form; 

(2it)  °°  J"  ’  •  •  •)  ^>1(k1-^1,k2-t2," '  ‘  )dB1d^2‘  •  • 

=  Jo(kltk2)  0  (rt)exp(-jk1  *r£)  +  (27r)  “  J-  •  -Jj  x(.^t  ,Bq  ,%2  »^2) 

C*>2*kl”’5o,^o-^l,k2,_^2’ "  * "  }  *  J(rt“Art  -  O  dAr^^d^d^'  •  ’ 

(28) 

The  til  ie  over  J  denotes  the  stochastic  Fourier  transform  (from  space 

to  ^  spacu),  J  is  the  transform  of  the  stochastic  part  of  the  physical 
no 

A  ->■  = 

optics  current  due  to  the  incident  field,  i.e.  2n  * H.  ,  and  x  is  3  dyadic 
kernel  involving  stochastic  transforms  of  Green’s  function  derivatives.  When 
the  illumination  is  unbounded  so  that  8  (r^)  =1  for  all  r^  ,  it  is  obvious 
that  (9)  is  valid.  However,  when  the  support  of  the  incident  illumination  is 
finite,  the  dependence  of  J  on  r^  is  very  difficult  to  obtain.  The  Fourier 
transform  of  (28)  with  respect  to  r  yields 


J"  j  J(i't.6i,l2,--*)<j51(k1-31,k242,---)dB1dI2---  =  J^.k^e^-k.  ) 

J  "  J  | So52  ^2  >  *2  (kl-60 '  V61  •  '  dird8o  | 


•  J(kt,ei,'52,*--)d61dl2--*  (29) 

■f 

It  is  tempting  to  assume  that  J  can  be  expressed  as  a  product  of  a  ^-dependent 
function  and  one  which  depends  upon  ’  However,  there  is  no  reason 

for  this  form  to  be  valid  in  general  because  Xfy  cannot  similarly  be  expressed 
as  a  product.  In  fact,  it  is  the  presence  of  the  integral  term  on  the  rhs  of 


(29)  which  greatly  complicates  matters.  This  term  is  important  when  multiple 
scattering  on  the  surface  is  significant;  in  the  absence  of  multiple  scatter- 
ing  J  has  the  same  r^-dependence  as  the  incident  illumination.  When  there 
is  strong  multiple  scattering  present,  the  ^-dependence  of  Jt  in  (29)  is 
influenced  by  the  statistical  properties  of  the  surface  through  cj> in  (29). 
This  result  has  a  simple  physical  explanation  in  that  multiple  scattering 
implies  that  the  current  at  any  point  on  the  surface  is  determined  by  the  behav¬ 
ior  of  the  surface  in  a  neighborhood  of  the  point;  thus,  the  support  of  the 
current  on  the  surface  will  differ  from  the  projected  area  of  the  incident 
illumination.  One  may  also  visualize  the  multiple  scattering  process  in  terms 
of  rays  undergoing  at  least  two  reflections  on  the  surface  before  traveling 
off  into  free  space;  consequently,  there  will  be  some  rays  reflected  to  the  sur¬ 
face  outside  of  the  incident  illumination  support.  What  this  all  means  is 
that  the  ^-dependence  or  the  support  of  the  transformed  current  is  determined 
both  by  the  incidence  illumination  and  the  surface  roughness  when  multiple 
scattering  is  important. 

To  overcome  this  difficulty  and  yet  not  reduce  the  problem  to  the  trivial 

case  of  single  scattering,  it  is  necessary  to  assume  that  the  support  of  the 

incident  illumination  projected  on  the  z*G  plane  is  very  large  compared  to 

A^  .  If  this  is  true  then  0(kt  -k^  )  will  be  very  sharply  peaked  about 

kt«k^  and  J  must  then  necessarily  have  the  same  dependence.  Physically, 

this  approximation  implies  that  the  illuminated  area  is  so  large  that  there  is 

little  relative  difference  between  the  projected  area  of  the  incident  wave  and 

the  true  illuminated  area  including  multiple  scattering.  Thus,  as  long  as  the 

cross  sectional  area  of  the  incident  beam  plane  wave  projected  onto  the  z  *  0 

2 

plane  is  large  compared  to  Aq  ,  (29)  shows  that 
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(30) 


J(rt,ki»K)  -  j(ltlfK)  6(rt)exp(-jk1  *rt> 

is  a  valid  approximation. 

Equation  (30)  is  a  very  important  result  but  one  must  be  careful  to  under¬ 
stand  what  it  does  and  does  not  imply.  Equation  (30)  can  be  used  to  compute 
the  statistical  moments  of  the  scattered  field;  it  cannot  be  used  to  infer  the 
sample-by- sample  behavior  of  the  scattered  field.  Thus,  just  because  J  has 
a  finite  support  does  not  mean  that  one  may  use  (2),  i.e.  the  far  field  approxi¬ 
mation,  rather  than  (1),  i.e.  the  exact  expression  for  ,  because  the  actual 
support  of  the  saiq>le-by-sample  current  may  vary  significantly  from  0(r^)  . 

This  is  a  point  that  must  be  remembered  when  attempting  to  make  far  field  mea¬ 
surements  of  the  field  scattered  by  a  rough  surface.  That  is,  if  one  bases 
the  far  field  condition  on  0(rt)  ,  it  must  be  realized  that  there  may  well  be 

-f 

some  samples  of  E  which  do  not  satisfy  this  condition.  Thus,  one  would  be 
s 

-f 

well  advised  in  a  measurement  program  to  exceed  any  distances  based  on  G(rt). 

As  it  turns  out,  the  statistical  moments  of  the  scattered  field  are  the 
same  regardless  of  whether  one  uses  (1)  or  (2)  in  conjunction  with  (30).  How¬ 
ever,  as  noted  above  there  is  no  rigorous  justification  for  starting  with  (2). 
Furthermore,  there  is  a  very  important  point  in  the  development  starting  with 
(1)  which  shows  why  the  far  field  approximation  fails  when  the  illumination  is 
unbounded.  To  simplify  matters,  the  development  will  be  carried  thru  for  the 
mean  scattered  field,  however,  the  principles  apply  also  to  all  the  moments  of 
the  scattered  field. 

Using  J(r)dr  for  3  (r)dS,  the  average  scattered  field  from  (1)  is 
c  s 

given  by 

<Es($o)>  “  *  J  ^  Vq  xVq  x  J < J (r ) G ( | RQ-r  |  ) >  d^  (31) 
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where 


;J(r)C(|Ro-?|)>  -J*  J  J(?t,C,a?)G(|Ro  -?t|  .  Zo-^)Pl(C,30dCdOO 


(32) 


Converting  the  integrals  in  (32)  to  convolutions  in  stochastic  Fourier  trans¬ 
form  space  yields  the  following; 

~(s„+l) 

JJJJ  I(?t,k1-60,I)G(|Rt 


<•>  -  (2it) 


o 


•  exp(jeoC)4)1(-k1,-K)dCdBodk1dK 


Substituting  from  (30)  yields; 


<•> 


(2u)  “  ^J'J'J'jCk^.K)  <X>1(-k1,-K)  6(rt)G(|Ro  ZQ-C) 


•  exp (-  j  k  *r  +  j3  OdCd6dk.dK 
t  o  or 


(33) 


Inserting  this  expression  in  (31)  and  rearranging  integrations  yields  the 
following  r^-integral; 

I  0(rt)G(|Ro  -rj  ,Zq-  C)exp ("j^  *rt)drt 
J  t  t 

which  may  be  rewritten  as  a  convolution  of  transforms,  i.e. 

(2tt)'2  I  e(-k±  +kt)G(kt  ,  Zo-r,)exp(-jkt4o  )dkt 
J  t  t 

-  (2ir)'2  J  0(kt-k^)e>T  -  WkZ-kf  |!0-{|](-J2VkJ-Jj)1^t 


(34) 


Taking  Zq  to  always  be  above  the  surface  Implies 


lZo  -  «l"Zo  -  «  • 


The  resulting  dependence  of  the  rhs  of  (34)  on  Rq  corresponds  to  a  generalized 
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plane  wave,  consequently,  the  curl  curl  operation  in  (31)  may  be  simplified 
to  the  following; 


VoxVox  (1 )  -  -  kQ2  1(1  (k  x  k  x  j) 


where  «  exp(- j  kt*RQ  -  jVk*-k*  Zq  ) /(-j2  )  and 


(k  +  Vk2  -  k*  z)  /  |k  +  Vk2  -  k2  z | 

'  t  o  t  /  /  1  t  o  t  1 


Substituting  these  results  in  (33)  and  (31)  yields 
k 


<E  (R  )> 
s  o 


f Uf  ^  **  X^(kl"6o’^  O  ) 


-1/2 


•  dktdUBodkdK/(2TT) 


s  +3 


(35) 


Since  for  any  realization  of  the  surface  Zq  >  C  *  the  dominant  terms  in  the 
^-integration  are  0  and  the  exponential  factor  containing  Rq  .  For  any 
9(rt)  having  finite  support  it  is  always  possible  to  take  Rq  large  enough 
that  the  k^-integral  can  be  asymptotically  evaluated  via  the  method  of  steepest 
descent  [Clemmow,  1966].  However,  if  0(r^)  has  infinite  support  then  it 
dominates  the  ^-integration  because  80^-iT^  )  =  6  0^-icj  )  .  Since  use  of 
the  steepest  descent  method  corresponds  to  the  far  field  approximation 
[Clemmow,  1966],  it  is  clear  that  the  validity  of  the  far  field  approximation 
is  entirely  dependent  upon  which  of  these  factors  is  dominant  in  the  £  -Integral. 
It  is  also  Interesting  to  note  that  the  k^-integration  is  essentially  of  the 
same  form  that  Hansen  and  Bailin  [1959]  studied  for  the  case  of  a  uniformly 


illuminated  deterministic  circular  aperture.  They  found  that  the  far  field 


approximation  is  valid  for  Rq  >  D  /Xq  where  D  is  the  diameter  of  the  aper- 


2  .  9  2  -*• 

ture.  Thus,  for  R  >  D  /A  and  ir(D/2)  »  A  ,  the  k  -integration  in  (35) 

o  ~  o  o  f." 


may  be  asymptotically  evaluated  using  the  method  of  steepest  descent.  Follow¬ 


ing  this,  one  completes  the  £  and  integrations  in  (35)  to  yield 


A  A 


<E  (R  )>  ‘  jk  n  G  (R  )0(k  -k  )k  xk  x 
so  oo  o  s^i^ss 


J J j(k1+kg  ,K)  <t1(-k1,-K)dk1dic/(2Tr)  ‘ 


(36 


where  k  is  a  unit  vector  in  the  direction  of  the  point  of  observation 


(R  ■  R  k  )  and  k  =  k  z  *  k  is  the  product  of  k  and  the  z-component  of 
oos  sos  r  o 

z 


this  unit  vector. 

This  same  asymptotic  approach  may  be  used  to  obtain  an  approximate  expres 
sion  for  the  variance  of  the  scattered.  To  account  for  the  limited  extent  of 


the  Incident  illumination,  it  is  sufficient  to  replace  j  by  j  0(r  )  and 

-*•-**->  t]L 
j*  by  j  0 (r  )  in  (18).  One  then  uses  the  steepest  descent  method  in 

2 


k-space  to  first  evaluate  the  r  -integral  and  then  the  r  -integral.  Asymp¬ 


totic  evaluation  of  the  r  -integral  leads  to  the  following  k-space  integra- 

C1 


tion; 


/ 


0(ks  -  kt  +k2)$2(-k2)exp(- j k2*rt  )dk2 


If  the  illuminated  area  encompasses  many  correlation  lengths  then  0  is  very 
sharply  peaked  relative  to  $>2  and  the  above  integral  may  be  approximated 


as  follows; 


*2(VV  J  ®(^st'itit+S)exp('jit2,?t2)di‘2  *  *2(VV 

0(rt  )exp[j(ic8  -ki  )  «rt  ] 
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Combining  Che  r  -dependent  terms  from  this  result  with  the  remaining  r 

2  C2 
factors  in  (18)  and  expressing  the  r  -integral  as  a  convolution  in  k-space 

C2 

yields  the  following  integral; 

02  (k  -k)G(-k,Z  _£)exp(jk*R  )dk 


fi 


When  this  integral  is  evaluated  by  the  steepest  descent  method,  one  obtains 
a  result  proportional  to  0^(0)  which,  for  uniform  illumination,  is  equal  to 
the  illuminated  area.  The  final  result  f_>r  the  variance  of  the  scattered  field 
is  identical  to  (16)  except  that  the  illuminated  area  A  is  finite.  Once 
again,  it  is  the  interplay  between  the  Green's  function  and  the  illumination 
support  that  determines  the  suitability  of  the  far  field  approximation. 

It  should  be  emphasized  that  these  results  are  very  general  in  that  they 

do  not  depend  on  any  special  approximation  for  the  current  induced  on  the  rough 

surface.  However,  these  results  are  only  valid  if  the  cross  sectional  area  of 

the  incident  illumination  beam,  as  projected  onto  the  z=0  plane,  is  very 

2 

large  compared  to  Aq  .  This  simplification  minimizes  the  error  due  to  the 
difference  between  the  illumination  area  and  the  true  area  on  the  surface  over 
which  the  current  is  nonzero.  Furthermore,  in  the  case  of  the  variance  of  the 
scattered  field,  it  is  also  necessary  to  stipulate  that  the  illuminated  area 
contain  very  many  decorrelation  intervals  of  the  surface  roughness  in  order  to 
recover  (16).  When  this  simplifying  requirement  is  not  met,  the  problem  becomes 
more  difficult  and  will  be  left  to  future  studies. 
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